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INTRODUCTION 


r 


Under  a previous  contract  ( F19F 28-73-C-O081)  with 
Emmanuel  College,  the  final  report  ( AFCRL-TR- 75-0588) 
has  detailed  the  specifications  of  the  mini-computer 
located  at  AFGL,  which  system  has  been  Installed  for 
data  reception  In  the  MAGAF  network. 

1 

Now  the  MAGAF  Data  Collection  system  consists  of  various 
programs,  functions,  and  subroutines  designed  to  collect, 
archive,  and  provide  user  access  to  data  received  from 
the  MAGAF  data  collection  network.  These  routines  operate 
under  the  VORTEX  real-time,  multltasked,  disk  operating 
system  on  a VARIAN  V-72  minicomputer.  (Multitasking  is 
a scheme  whereby  many  more  or  less  Independent  programs 
compete  for  system  resources  under  some  kind  of  priority 
setup.)  A brief  description  of  MAGAF  system  flow 
fol 1 ows . 

Upon  system  Initialization,  a MAGAF  Initialization 
program  Is  read  In  from  the  Foreground  Library  and  executed. 

Since  this  program  Is  on  disk.  It  can  be  changed  as  required. 

Currently,  this  program  checks  the  presence  of  the  5 kHz 
and  10  kHz  signals  to  V-72  counter  and  Interrupt  clock  and 
syncs  the  computer  and  Its  counter  to  the  10  second  pulse 
from  the  master  clock.  At  this  point,  control  Is  returned 
to  the  VORTEX  operating  system  with  the  core  resident  DDU 
monitor  active. 

4 
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The  core  resident  monitor  initiates  a read  into  the 
first  of  two  buffer  areas  and  releases  control  to  the  VORTEX 
system.  Upon  receipt  of  a frame  of  data  from  the  DDU, 

VORTEX  reactivates  the  monitor  which  updates  the  pointer 
to  the  data  for  user  programs  and  initiates  a read  into  the 
second  buffer  area.  This  procedure  continues  using  alternate 
buffer  areas.  Subroutines  are  available  on  the  system  library 

to  access  the  data,  convert  the  time  code,  and  unpack  the 
data  into  Instrument  order. 

The  data  access  monitor  subroutine,  when  called  by 
a user  program,  checks  for  a change  in  the  data  pointer 
provided  by  the  DDU  monitor.  If  there  has  been  no  change 
in  this  pointer,  the  program  is  deactivated  for  50  milli- 
seconds and  the  pointer  Is  checked  again.  If  the  pointer 
has  changed,  the  frame  is  moved  to  the  data  area  defined 
by  the  user  program  and  the  subroutine  returns  to  the  user 
program.  It  should  be  noted  at  this  point  that  a change 
In  this  procedure  Is  contemplated. 

The  Interim  archiving  procedure  takes  the  frame  of  data, 
compresses  It  when  possible,  stores  It  on  disk  In  a circular 

buffer  to  maintain  a short  term  history,  and  schedules  the 

0 

tape  writing  program  when  a buffer  segment  has  been  completed. 
Both  programs,  compression  and  tape  writing  are  on  disk  and 
can  therefore  be  modified. 

In  the  future,  it  is  expected  that  among  other  thlnos 
the  following  items  will  be  implemented: 

Software  lock  on  archive  tape  and  protect  of  tape 
Software  system  reset  to  protect  archiving 
Network  reliability  on  a by-station  basis 
Graphics  capability  - real  time  and  from  tape 
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1.  MAGAF  DATA  PROCESSING  SYSTEM  ADDITIONS 


(a)  The  library  functions  listed  in  Appendix  A have  been  added 
to  the  system. 


Function 

Name 

Type 

Cal  1 ing 
Sequence 

Resul t 

IAND 

Integer 

i*iand(j,k) 

Returns  logical  product  of  J and  K. 

I0R 

II 

I*I0R(J,K) 

Returns  logical  sum  of  J and  K. 

IXOR 

II 

1*1 X0R(J , K) 

Returns  exclusive  of  J and  K. 

NOT 

II 

I-NOT(J) 

Returns  logical  complement  of  J. 

1 BIT 

II 

i*ibit(j,n) 

I={?)  if  Nth  bit  (mod.  16)  of 
J is  {°}. 

MASK 

It 

l-MASK(N) 

S’ N leading  ones  (>^\ 

Re,urns{|»|  l«din9  fIf"|<o| 

C zeros_J  V.  J 

ISHIFT 

II 

l-ISHIFT 

(J.N.ITP) 

Returns  J shifted  N places  according 
to  ITP. 

0 - Arithmetic  left  shift 

1 - Logical  rotation  left 

2 - Aritmetic  right  shift 

3 - Logical  right  shift 

(0  fill) 

1 SET 

II 

l-ISET(N) 

Returns  a mask  with  bit  N set  to  1 
(others  0).  N out  of  range  0-15 
undefined. 

IRSET 

II 

l-IRSET(N) 

Returns  a mask  with  bit  N set  to  0 . 

(others  1).  N out  of  range  0-15  is 
undefined  except, 

N = 17  returns  left  byte  mask 
N * 1 8 " right  byte  mask 

Li 





MAGAF  data  collection  interim  archiving  procedure 


The  MAGAF  interim  archiving  procedure  consists  of  two  programs  which  prepare 
the  archive  tape  and  a library  subroutine  and  program  for  using  previously 
prepared  tapes.  A description  of  the  programs  follows.  Program  listings 
are  contained  in  Appendix  B.  The  tape  format  is  given  in  Appendix  C. 


I.  Program  ARCH IV  monitors  the  incoming  data,  compresses  it  if  possible, 
and  buffers  it  out  to  disk.  When  the  disk  buffer  is  full,  ARCHIV  schedules 
program  ARCTAP,  which  transfers  the  data  to  tape  and  changes  tape  units 
when  the  current  unit  is  full. 


II.  Function  l$TPY  returns  tape  status  and  frame  data  to  user  programs, 
scheduling  program  DCMPRS  when  necessary.  Program  DCMPRS  reads  in  one 
physical  record  from  tape,  decompresses  the  data,  writes  it  out  to  disk 
and  returns  to  l$TPY  tape  status  and  the  number  of  logical  records 
contained  to  l$TPY. 

The  calling  sequence  for  l$TPY  is: 


where  1ST  is  tape  status,  IUNIT  is  the  logical  unit  number  of  the  tape 
and  IA  is  an  array  with  dimension  2A5  which  contains  frame  and  status 
information. 

1 - good  frame 
0 - device  busy 
1ST  = -1  - end  of  tape 

-2  - end  of  file 
-3  " error 


The  format  of  IA  is  as  described  in  AFCRL-TR-75-0588 , with  the  following 
exceptions : 


A.  Data  words  and  buffer  contents  are  in  the  following  form 


sign  bit  (repl icated) 
data  bit 
error  bit 


where  s 


B.  IA(242)  and  I A (243)  now  contain  frame  time  in  seconds  and  mi  1 1 i 
seconds  respectively,  frame  time  being  reset  to  zero  every  time 
real  time  reaches  zero  modulo  ten  seconds. 
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000002  P CALL 


COMPLETE 


000000  R ISHIFT  ©00000  R RETU 


VORTEX  DASMP 


eeis  hours 


1 

2 

3 


4 

NAME 

ISET 

O0C00O 

ooooeo 

A 

5 

I SET 

ENTP 

6 

* 

7 

» 

3 

T 

000000 

Q 

RETU 

BES 

0 

10 

EXT 

*SE 

000001 

002000 

A 

1 1 

CALL 

*S£ 

000002 

00O0OE 

000000 

E 

12 

CALL 

BES 

0 

OOO0O3 

000001 

A 

13 

DATA 

1 ,0 

000004 

eooooo 

A 

000005 

074004 

A 

14 

STX 

5TX+1 

000006 

0C6037 

A 

IS 

LDXE* 

CALL +2 

000007 

100064 

R 

000010 

015421 

A 

16 

LDA 

0421 , 1 

000011 

00603O 

A 

17 

STX 

LDXI 

* 

000012 

00001 1 

P 

0.00013 

001000 

A 

18 

RETU* 

RETU 

C00014 

100000 

R 

19 

END 

ENTRV  NAMES 
OOOOOO  P I SET 
EXTEPNAL  NAMES 

oooooe  f tsr 

SYMBOLS 

000002  E »SE  000002  R CALL  000000  P 
000011  R STX 

0 ERRORS  ASSEMBLY  COMPLETE 


J* ISET ( N) 

RETURNS  A EXT  MASK.  IJITH  PIT  N •= 
1 (OTHERS  O'.  N OUT  OF  RANGE  0- 


ISET  000000  R RETU 


VORTEX  DASMR 


001S  HOURS 


1 

NAME 

IRSET 

000000 

000000 

A 

2 IRSET 

3 * 

4 * 

5 * 

6 * 

SMTP 

000000 

7 RETU 

BES 

0 

3 

EXT 

*SE 

000001 

002000 

A 

9 

CALL 

*SE 

000002 

000002 

000000 

E 

10  CALL 

BES 

0 

O0O003 

000001 

A 

11 

DATA 

1.0 

000004 

ooooeo 

A 

000005 

074004 

A 

12 

STX 

STX+1 

000006 

006037 

A 

13 

LDXE* 

CALL+2 

000007 

100004 

R 

000010 

015441 

A 

14 

LDA 

0441 , 1 

000011 

0O603O 

A 

15  STX 

LDXI 

* 

000012 

000011 

R 

000013 

O010OO 

A 

16 

RETU* 

RETU 

000O14 

100000 

R 

17 

END 

ENTRV  NAMES 


000000  R 
EXTERNAL 
000002  F 
SYMBOLS 
OOOOOr.  F. 
0000  I I R 


IPSET 

NAMES 

*SE 

»3F 

STX. 


J-IRSET(N) 

RETURNS  A BIT  MASK  LIITH  BIT  N SE 
0 (OTHEPS  1 > . N OUT  OF  RANGE  0-1 
EXCEPT  17  GIVES  LEFT  BYTE  MA3F  . 1 


AOOO02  R CALL  000000  R IPSET  000000  R RETU 


O f PROPS  ASSEMBLY  COMPLETE 


VORTEX  DASMR 


0000  HOURS 


000422 

A 

1 

TWO 

eou 

0422 

0O0423 

A 

2 

FOUR 

ECU 

0423 

000424 

A 

3 

EIGHT 

EOU 

0424 

000466 

A 

4 

SIX 

EOU 

0466 

O0046? 

A 

5 

SEVEN 

EOU 

0467 

£100443 

A 

6 

115 

EOU 

0443  MINUS  FIVE 

000465 

A 

7 

FIVE 

EOU 

0465 

000464 

A 

3 

THREE 

EOU 

0464 

000472 

A 

9 

FIFTH 

EOU 

0472 

000002 

A 

10 

B 

EOU 

2 

000001 

A 

11 

x 

EOU 

1 

000025 

A 

12 

TUNIT 

EOU 

21 

000041 

A 

13 

DUNIT 

EOU 

33 

14 

ARC 

REU 

FCB, DUNIT 

000000 

006505 

A 

000OO1 

OO0O0O 

E 

• 

00OO05 

100000 

A 

000003 

001441 

A 

000004 

000175 

R 

00OO0S 

0OO0OO 

A 

O0OOO6 

000000 

A 

00000? 

006010 

A 

15 

LDAI 

TUNIT 

000010 

000025 

A 

000011 

050005 

A 

16 

STA 

S 

00001c.' 

006010 

A 

17 

LDAI 

DUNIT 

000013 

000041 

A 

000014 

050006 

A 

18 

STA 

6 

000016 

0060 1 0 

A 

19 

LDAI 

FCB 

000016 

000175 

P 

O0OO1? 

050004 

A 

20 

STA 

4 

000050 

OO2O0O 

A 

21 

ARCH1 

CALL 

ARCHIV 

000021 

0OO211 

R 

O0OO22 

12414? 

A 

22 

ADD 

NUD1 

00O023 

005012 

A 

23 

TAP 

000024 

144146 

A 

24 

SUB 

DCF 

000026 

001004 

A 

25 

JAN 

ADD 

O0OO26 

000067 

P 

00002? 

005301 

A 

26 

DECP 

01 

000030 

002000 

A 

27 

CALL 

OUT 

00OO31 

000075 

R 

000032 

014145 

A 

23 

LDA  • 

FCB  + 3 

O0OO33 

124154 

A 

29 

ADD 

TUOTUO 

000034 

144152 

A 

30 

SUB 

NSECT 

O0OO35 

054142 

A 

31 

STA 

FCB  + 3 

O0O03 C 

005001 

A 

32 

TZA 

00003? 

05414? 

A 

33 

STA 

NSECT 

000040 

014131 

A 

34 

LDA 

NIJD1 

35 

SCHED 

31 , 1 , 106 , 'F ' , ' AR ' , ' CT ' , ' AP ' 

00004 1 

006505 

A 

000042 

OOOOO0 

C 

000043 

01013? 

A 

00OO44 

145152 

A 

000046 

140722 

A 

VORTEX  DASMR 


0000  HOURS 


■ 


! 


i 


5 

4 

H 

h 


000046 

141724 

A 

00004V 

140720 

A 

000050 

024524 

A 

36 

LDB 

MUDS 

00005 L 

014126. 

A 

37 

LDA 

FCB  + 3 

000052 

12  1135 

A 

33 

ADO 

TUOTUO 

000053 

144125 

A 

39 

SUB 

FCB+4 

O0OO54 

001004 

A 

40 

JAN 

ADD 

000055 

000067 

p 

00OO56 

001010 

A 

41 

JhZ 

ADD 

00005? 

000067 

R 

42 

REUI 

FCB. DUNIT 

000060 

005505 

A 

00006 1 

OOOOOl 

E 

000062 

100000 

A 

00OO63 

001441 

A 

00O064 

0001 75 

P 

000065 

000000 

A 

000066 

000000 

A 

000067 

064102 

A 

43  ADD 

STE 

NUD1 

000070 

014504 

A 

44 

LDA 

NWDS 

000O71 

0020O0 

A 

45 

CALL 

OUT 

000072 

OO0075 

R 

000073 

0O10OO 

A 

46 

JMP 

ARCH1 

00OO74 

000020 

R 

47  * 
43  X 
«*0  * 

000075 

000OC0 

A 

50  OUT 

ENTP 

000076 

054464 

A 

51 

STA 

TUP 

000077 

006057 

A 

52 

STAE* 

ABUF 

000100 

100171 

R 

000101 

001004 

A 

53 

JAN 

WRITE 

000102 

0001 15 

R 

000103 

014453 

A 

54 

LDA 

HDAD 

000104 

054454 

A 

55 

STA 

STRT1 

OO01OS 

01 4061 

A 

56  LPT  1 

LDA 

LEFT 

O00106 

00501 4 

A 

S7 

TAX 

00010V 

144453 

A 

58 

SUB 

TMP 

000110 

00531 1 

A 

59 

OAR 

0001 1 1 

00 1 002 

A 

60 

JAR 

FILL 

O00112 

000144 

R 

000113 

006506 

A 

61 

JSR 

FILLIN, B 

000114 

00O156 

R 

62  WRITE 

WRITE 

FCB, DUNIT 

000115 

OO6505 

A 

000 1 1 6 

000061 

E 

0001 IV 

100000 

A 

O0O120 

00O44J 

A 

000121 

000 1 75 

R 

000122 

OOOO00 

A 

00O1  2 -i 

000000 

A 

000121 

04  H'f.2 

A 

r.3 

1 NR 

N5ECT 

0001  ftU 

014  4 < : 

A 

6 1 

t.l'A 

5TWT1 

18 


000126 

1.  4040 

0 

06 

A PI) 

If  r T 

OOu  1 1; 

0614  : 1 

A 

t . 

•.'I'O 

•.TRT  1 

OP o 1 ";o 

oi  i4:.',: 

.3 

i.7 

L I)i  i 

TMP 

oeioi'Ji 

1 4 40  I1. 

1 1 

iy.': 

611 1; 

LF.  i T 

eooi 3a 

0O4O4 1 

A 

' .9 

L.  0 1-: 

DC  Hi  1 

000133 

064035 

A 

70 

STB 

ABUF 

000134 

024033 

A 

71 

LDB 

D120 

000135 

064031 

A 

72 

STB 

LEFT 

000136 

054424 

A 

73 

STft 

TMP 

000137 

005311 

A 

74 

DAP 

000140 

001002 

A 

75 

JAP 

LFT 1 

000141 

000105 

R 

OO0142 

001000 

A 

76 

RETU* 

OUT 

O0O143 

100075 

R 

000144 

005111 

A 

77  FILL 

IAR 

000145 

054021 

A 

7S 

STA 

LEFT 

000146 

034414 

A 

79 

LDX 

TMP 

000147 

006506 

A 

SO 

JSR 

FILLIN, B 

000150 

000156 

R 

000151 

014411 

A 

31 

LDA 

TMP 

000152 

124016 

A 

82 

ADD 

ABUF 

000153 

054015 

A 

33 

STA 

ABUF 

000154 

001 000 

A 

84 

RETU* 

OUT 

000155 

100075 

R 

35  X 

36  X 

77  X 

000156 

005344 

A 

??■  FILL  IN 

DXP 

000157 

006215 

A 

39 

LDAE* 

STRT1 , X , 0200 

000160 

100561 

R 

000161 

006255 

A 

90 

STAEX 

ABUF, X, 0200 

000162 

100171 

F: 

000163 

001046 

A 

91 

JXNZ 

FILLIN 

000164 

000156 

R 

000165 

00S706 

A 

92 

IJMP 

O , B 

000166 

oooooo 

A 

93  * 

94  *. 

95  X 

0001C7 

000170 

A 

96  LEFT 

DATA 

120 

000 170 

000170 

A 

97  D120 

DATA 

120 

000171 

001 253 

R 

92  ABUF 

DATA 

BUFR 

000172 

oooooo 

A 

99  nuni 

DATA 

0 

100  DCE 

DIP 

2660 , BUF  R 

0001 73 

00'  .01 .10 

A 

0 0 l 74 

00 1 253 

R 

OOO 1 75 

oooooo 

A 

000175 

lot  fob 

OPG 

*-l 

102 

FOB 

120, BUFR, 3 

000 1 75 

OOO 1 70 

A 

OOO  1 7 Ci 

001263 

P 

000177 

0 0 1400 

A 

O0O20O 

OOOOOO 

A 

VORTEX  DASMR 


9009  HOURS 


O00S01 

000000 

A 

000302 

000000 

A 

000203 

000000 

A 

0OO204 

000000 

ft 

O00205 

000000 

A 

OO02OS 

000000 

A 

000207 

000000 

A 

103 

NSEC  T 

DATA 

0 

000210 

000026 

A 

104 

TUOTUO 

DATA 

22 

105 

X 

106 

* 

107 

* 

000211 

R 

10S 

arch:v 

EOU 

* 

000211 

000212 

R 

109 

DATA 

*+l 

110 

DELAY 

DELAY 

1 

00O212 

0O6505 

A 

000213 

000042 

E 

000214 

001100 

A 

00O215 

000001 

A 

000216 

000000 

A 

000217 

010002 

A 

111 

DEI.  VI 

LDA 

2 

000220 

144351 

A 

112 

SUB 

CUR 

O00221 

001010 

A 

113 

JAZ 

DELAY 

00O222 

000212 

R 

O00223 

030002 

A 

114 

LDX 

2 

000224 

074345 

A 

115 

STX 

CUR 

O00225 

015364 

A 

116 

LDA 

244.X 

0OO226 

054346 

A 

1*7 

STA 

NUDS 

O0O227 

015360 

H 

113 

LDA 

240.X 

O00230 

054345 

A 

119 

STA 

STAT 

000231 

015363 

A 

120 

LDA 

243.X 

000232 

054356 

A 

121 

STA 

BUFF 

000233 

02OOO3 

A 

122 

LDB 

3 

TIME 

000234 

005001 

A 

123 

TZA 

090235 

004441 

A 

124 

LLRL 

' 1 

0O0236 

004141 

A 

125 

LSR8 

1 

000237 

006170 

A 

126 

DIVX 

5000 

000240 

011610 

A 

000241 

064345 

A 

127 

STB 

SEC 

SECONDS 

O0O242 

120464 

A 

128 

ADD 

THREE 

ROUND 

000243 

00456O 

A 

129 

LLSR 

16 

0OO244 

170465 

A 

130 

DIV 

FIVE 

000245 

064342 

A 

131 

STB 

MILLI 

MILLISECONPS 

132 

EXT 

•UNPK 

000246 

002000 

A 

133 

CALL 

• UNPK , (CUR)*. 

sex 

000247 

000000 

E 

00O250 

100572 

P 

090251 

000012 

R 

000252 

034322 

A 

134 

LDX 

NUDS 

000253 

0051 44 

ft 

l 35 

IXR 

OO0254 

00534 4 

A 

136 

DXR1 

D.VR 

009255 

00f.O  l «; 

ft 

137 

LDftE 

BUFF , X 

0002'...:. 

000611 

ft 

000257 

004501 

A 

13C 

LflSP 

1 

20 


I 


■ 


\ „ ’ 

VORTEX  DASMP  0000  HOURS 


OOCSGO 

004304 

A 

139 

ASRA 

4 

OO02G1 

004401 

A 

140 

LASL 

1 

000363 

006065 

A 

141 

STAE 

BUFF  , X 

0OO363 

0006 1 1 

R 

000364 

001046 

A 

142 

jxnz 

DXR1 

O0O365 

000254 

R 

000366 

014307 

A 

143 

LDA 

STAT 

000367 

001002 

A 

144 

JAP 

PACK 

000370 

000307 

R 

000371 

0O6O1O 

A 

145 

LDAI 

SEC 

000373 

0 0 0 6 0 7 

P 

C0O273 

054265 

A 

146 

STA 

STRT1 

000274 

006010 

A 

147 

LDAI 

IPACK 

000275 

000577 

P 

000276 

054261 

A 

143 

STA 

START 

000277 

010464 

A 

149 

LDA 

THPEE 

000300 

124274 

A 

ISO 

ADD 

MUDS 

000301 

054002 

A 

151 

STA 

* + 3 

000302 

006505 

A 

152 

JSR 

MOV,  X 

0OO3O3 

000540 

P 

000304 

OO50O0 

A 

153 

NOP 

O003O5 

0O1O0O 

A 

154 

JMP 

RET1 

000306 

000353 

P 

000307 

P 

155  PACK 

EOU 

* 

156 

EXT 

SNAP 

000307 

006010 

A 

157 

LDAI 

sex 

00031 0 

OOOC 1 2 

P 

0 

000311 

054246 

A 

153 

STA 

START 

000312 

054246 

A 

153 

STA 

SI  PT 1 

000313 

006010 

A 

160 

LDAI 

IPK1 

000314 

000577 

R 

000315 

054244 

A 

161 

STA 

I PKURT 

000316 

006010 

A 

162 

LDAI 

. 50 

000317 

O0O062 

A 

000320 

002000 

A 

163 

CALL 

CMPR5 

000321 

0007)60 

R 

000322 

0OS0OO 

A 

164 

CALI. 

CMPRS 

000323 

000 '.60 

R 

000324 

002000 

A 

165 

CALI. 

CMPRS 

000  325 

0007160 

R 

000326 

0060 1 0 

A 

166 

LDAI 

10 

O00327 

0000  I 2 

A 

0003 30 

002000 

A 

167 

CALL 

CMF'RS 

000331 

000  36.O 

P 

000332 

06*2000 

A 

163 

CALI. 

CMPRS 

0O‘>  3 33 

000 360 

P 

0007)34 

06*2000 

A 

169 

CALL 

CMPRS 

000335 

000360 

P 

000336 

006505 

A 

170 

JSR 

MOV,  X 

000337 

000540 

R 

000340 

00001  1 

A 

171 

DATA 

9 

000341 

006010 

A 

172 

LDAI 

23 

000342 

000027 

A 

21 


W&r  * 


! 

t 


' 


VORTEX  DASMK* 


000343 

005000 

A 

173 

CALL 

CMPRS 

000344 

000360 

R 

O0O345 

010472 

A 

174 

LDA 

FIFTN 

000 346 

002000 

ft 

175 

CALL 

CMPRS 

00O347 

00O360 

9 

000350 

006505 

A 

176 

JSP 

MCV,X 

000351 

00054O 

ft 

000352 

00O015 

ft 

177 

DATA 

13 

000353 

014204 

A 

173 

RET1 

LDA 

START 

000354 

144202 

A 

179 

SUB 

HDAD 

000355 

054217 

A 

120 

sta 

NUDS 

000356 

001000 

A 

121 

RETU7 

ARCHIV 

000357 

10021 1 

P 

122 

t 

123 

* 

124 

% 

000360 

000000 

A 

185 

C-MPRS 

EHTR 

000361 

005311 

A 

186 

DAP 

000362 

054144 

A 

127 

STA 

N 

000363 

006027 

rt 

1 85 

LPBF* 

STPT1 

000364 

100561 

R 

000365 

002000 

A 

129 

CALL 

INSRT 

O0036S 

000532 

R 

000367 

005311 

A 

190 

DAP 

000370 

006O14 

A 

191 

TAX 

O0O37J 

124167 

A 

1 ‘ ? 

ADD 

STRT1 

000372 

005012 

A 

i<: 

TAB 

O00373 

006001 

A 

194 

T2A 

000374 

064166 

A 

195 

STA 

TMP 

000375 

016001 

A 

i 96 

LDA 

LDA 

1 . B 

000376 

146000 

A 

197 

SUB 

0,  B 

O0O377 

004301 

A 

198 

ASPA 

1 

000400 

OO6055 

A 

199 

STAE 

- TEMP. X 

000401 

001172 

R 

000402 

001002 

A 

200 

JAP 

SUB 

000403 

000410 

P 

000404 

005211 

A 

201 

CPA 

000405 

001015 

A 

202 

JANZ 

* + 3 

O0O4O6 

000410 

R 

000407 

005111 

A 

203 

IAR 

000410 

144152 

A 

204 

SUB 

SUB 

TMP 

000411 

001004 

A 

205 

JAN 

*44 

000412 

000415 

P 

O0O413 

124147 

A 

206 

ADD 

TMP 

00041 4 

054 1 46 

A 

207 

STA 

TMP 

000  41 r. 

00". .322 

ft 

208 

DPR 

00041 6 

006  3-14 

ft 

809 

DXR 

000417 

001046 

A 

2 1 0 

JXNZ 

l-PA 

000420 

000  "J  75 

h’ 

ooo-i, :i 

010443 

ft 

8t  i 

l.PA 

MS 

000422 

024 1 40 

A 

218 

LDP 

TMP 

000423 

On]  0,30 

A 

21  N 

JPZ 

DONE 

000424 

000441 

R 

22 


• % 
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VORTEX  DASMR  0000  HOURS 


000425 

005111 

A 

214 

IAR 

IAP 

000426 

001010 

A 

215 

J AZ 

DN1 

000427 

000435 

p 

0 0043  0 

001020 

M 

216 

JBZ 

DONE 

000431 

000441 

R 

000432 

004101 

M 

217 

ASRB 

1 

000433 

00 1 000 

fi 

21 S 

JMP 

IAR 

000434 

000425 

R 

O0O435 

00511 1 

A 

213 

DN1 

IAR 

000436 

004102 

A 

220 

ASPB 

2 

000437 

001Q26 

H 

221 

JBNZ 

TUEL 

O0O4  40 

000521 

R 

000441 

120465 

A 

222 

DONE 

ADD 

FIVE 

0 0 C 4 4 2 

006057 

A 

223 

STAE* 

IPKURT 

000443 

100562 

R 

000444 

054121 

A 

224 

STA 

MBITS 

000445 

OO101O 

H 

225 

JAZ 

ZEITS 

000446 

000510 

p 

000447 

00500 1 

f\ 

226 

TZA 

000450 

054114 

A 

227 

STA 

N1 

000451 

0241 16 

H 

223 

LDB 

TWELVE 

000452 

174113 

A 

229 

DIV 

NBITS 

000453 

064113 

A 

230 

STB 

NTOUD 

000454 

014114 

A 

231 

LDA 

LLSR 

000455 

124110 

A 

2 72 

ADD 

NBITS 

000456 

054007 

A 

STA 

LLS 

000457 

006010 

A 

2 ' 4 

LDAI 

TEMP 

O0O4CO 

0O 1 t 72 

R 

000461 

064)02 

fi 

235 

STA 

TMP I 

O00462 

034104 

H 

236 

LDX 

NTOUD 

000463 

0060 1 7 

H 

2 37 

NXT1 

LD«E* 

TMP1 

000464 

100564 

P 

000465 

044076 

A 

238 

I NR 

TMP  1 

000466 

004540 

H 

2 39 

r 

r 

cn 

LLSR 

0 

0OO467 

0O5344 

A 

240 

DXR 

00O470 

001046 

A 

24  1 

JXHZ 

NXT 1 

000471 

000463 

P 

000472 

004544 

A 

24  2 

LLSR 

4 

O0O473 

002000 

A 

243 

Call 

I NSRT 

00047 4 

00*1632 

R 

000475 

01406V 

H 

24  4 

LDA 

HI 

000476* 

12407 0 

A 

246 

ADD 

NTOUD 

000477 

064066 

A 

24fci 

STA 

HI 

000500 

1 44076 

A 

£47 

SUE 

N 

00050 1 

001004 

A 

£48 

JAN 

NXT 1-1 

000302 

000462 

P 

0OO5O3 

0714023 

H 

249 

LDX 

N 

0007*04 

006276 

A 

260 

LDBE* 

STRT1 . X , 0200 

000506 

1 00'.  *'.  1 

P 

0005 OC 

007.000 

A 

251 

CALL. 

1NSRT 

OOO'  .07 

000632 

f: 

000510 

O 1 4060 

M 

767 

xniri. 

Ll'A 

6TRT1 

0(105  1 1 

1.:  i0|6 

A 

26  < 

ADD 

N 

23 


VORTEX  PASMR  0000  HOURS 


000512 

005111 

0 

254 

IAR 

0OO513 

054045 

255 

sta 

STRTl 

00051  4 

014O12 

0 

256 

LODX 

LDA 

N 

000515 

005111 

A 

2P7 

IAP 

000516 

044043 

A 

253 

IMP 

IPKUPT 

000517 

001000 

A 

253 

RETU* 

ompps 

O0052O 

100360 

p 

000521 

034046 

A 

260 

TWEL 

LDX 

TWELVE 

O0O522 

006077 

A 

261 

STXE* 

IPKURT 

000523 

100562 

R 

000524 

044034 

A 

262 

INP 

STRTl 

O0O525 

006505 

A 

263 

JSR 

MOV , X 

000526 

000540 

P 

00O527 

000000 

A 

264 

N 

DATA 

0 

0OO530 

001000 

A 

265 

JMP 

LODX 

000531 

000514 

R 

000532 

OOOOO0 

A 

266 

INSRT 

ENTR 

000533 

006067 

A 

267 

STBE* 

START 

O0O534 

1OO560 

P 

000535 

044022 

A 

£65 

INR 

START 

O0O536 

001000 

A 

269 

RETU* 

INSPT 

000537 

1O0532 

R 

270 

t 

27 1 

¥ 

272 

* 

O0O540 

025O0O 

A 

2*^3 

MOV 

LOB 

0,X 

000541 

005144 

A 

’ * 

IXR 

00O542 

074013 

A 

2'  5 

STX 

JUMP+1 

000543 

03401 4 

A 

276 

LDX 

START 

000544 

005322 

A 

277 

DERI 

PER 

000545 

006017 

A 

278 

LPAEX 

STRTl 

000546 

100561 

R 

00054? 

04401  1 

A 

279 

INR 

STRTl 

00O55O 

055000 

A 

2S0 

STA 

0.X 

00O551 

005144 

A 

23 1 

IXR 

000552 

00 1026 

A 

2-22 

JBNZ 

DBP1 

000553 

0OO54-) 

P 

00O554 

074003 

A 

233 

STX 

START 

000555 

00 1 000 

A 

264 

JUMP 

JMP 

t 

000556 

000555 

P 

285 

* 

as  e 

237 

t 

000557 

000573 

R 

233 

HD  AD 

DATA 

HEAP 

000560 

00OO0O 

A 

289 

START 

DATA 

0 

000561 

000000 

A 

£90 

STRTl 

DATA 

0 

0005.6? 

000000 

A 

2'4 1 

IPKI.IRT 

DATA 

0 

000563 

000000 

1*1 

292 

TMF' 

D A T A 

0 

0005.6  4 

00 110 00 

A 

29  3 

TMP1 

DATA 

0 

000505, 

000001 1 

A 

294 

Ml 

DATA 

0 

00^' . i . r. 

01111111)0 

A 

295 

Njiirr: 

DATA 

0 

Hf)»V 

01)0000 

A 

•V’li. 

NT0i.ll' 

DATA 

0 

OOoV/O 

0000 | 4 

A 

nit  i vfc 

Data 

12 

. 


I 


24 


VORTEX  DASMR  0000  HOU‘ 


000571 

004540 

292 

LLSR 

LL5P 

0 

00057c 

00*0000 

A 

£39 

CUP 

DATA 

0 

300 

* 

301 

X 

30  £ 

* 

000573 

000001 

A 

303 

HEAP 

DATA 

1,0 

000574 

000000 

A 

000575 

000000 

A 

304 

NUDS 

DATA 

0 

000576 

305 

STAT 

BSS 

1 

000577 

000000 

A 

306 

I PACK 

DATA 

0 

000577 

R 

307 

IPK1 

EOU 

I PACK 

000600 

303 

BSS 

000607 

309 

SEC 

BSS 

1 

000610 

310 

MILL  I 

BSS 

1 

000611 

311 

EUF  r 

BSS 

1 

000612 

31£ 

SCX 

BSS 

£40 

001 172 

313 

TEMP 

BSS 

49 

001253 

314 

BUFR 

BSS 

ISO 

000000 

R 

315 

END 

ARC 

ENTRY  NAME'S 

external  manes 


000247 

E 

SIJNPK 

000000 

E 

SNAP 

000213 

E 

VSEXEC 

000116 

E 

V*IOC 

SYMBOLS 

000247 

E 

SUNPK 

000171 

R 

ABUF 

OO0067 

R 

ADD 

000000 

P. 

ARC 

j O0O02O 

P 

AF0H1 

OO021 1 

P 

ARCH IV 

000002 

A 

B 

00061 1 

R 

BUFF 

001253 

R 

BUFP 

000360 

p 

CMPP3 

00 OS 7 2 

F: 

CUR 

000170 

P 

0120 

000544 

P 

DBP 1 

OOO 173 

R 

DCB 

O0O2 l 2 

R 

DELAY 

000817 

R 

DELY1 

000435 

R 

PHI 

000441 

R 

DONE 

00004 1 

A 

DHNIT 

00025-1 

R 

L'XRl 

O0O424 

A 

FIOHT 

0001 75 

P 

FCfci 

000472 

A 

FIFTH 

000 144 

P 

FILL 

0C0156 

R 

F ILL  IN 

000465 

A 

FIVE 

O00483 

A 

F OUR 

O00567 

R 

HOAD 

000573 

P 

HEAD 

000425 

R 

I MR 

00 O 6 38 

R 

IIISPT 

0005V 7 

F' 

IF'ACK 

000577 

P 

IPK1 

000562 

P 

JPKURT 

0 0 0 5 «.  6 

P 

JUHF' 

OO0375 

LOrt 

000167 

P 

LEFT 

000106 

P 

LFT1 

00 04 06 

f: 

LI-5 

0OO571 

R 

LLSR 

O006 ! 4 

P 

LCD* 

000443 

A 

M5 

0006 1 0 

R 

MILL  I 

000540 

p 

MOV 

000587 

R 

N 

000565 

R 

N 1 

000566 

R 

MBITS 

000207 

P. 

NSECT 

' 000567 

R 

NTOUP 

00017S 

R 

NUD1 

O0057S 

R 

NUD5 

0G0463 

R 

NXT1 

0 0 & & 7 5 

P 

OUT 

O00307 

R 

PACK 

O 0 0 3 5 3 

R 

RET1 

©00612 

R 

sex 

A 0 0 G 0 7 

P 

SF  C 

000467 

A 

SEVEN 

000466 

A 

six 

000000 

E 

SNAP 

0 0 0 S '£>  ^ 

P 

START 

000576 

P 

STAT 

000661 

R 

STPTl 

0004 1 0 

R 

SUB 

001172 

P 

TEMP 

000464 

A 

THREE 

000563 

R 

TMP 

000564 

R 

TMP  1 

5 OOOtfeS 

A 

TUNIT 

$00521 

R 

TUEl. 

O0O57O 

P 

TWELVE 

O00422 

A 

TUO 

C0021 0 

R 

TWO  TIJO 

000213 

F 

VCEXEC 

000116 

E 

VHOC 

0001 15 

R 

WRITE 

0 0 0 0 0 1 

A 

*/ 

006510 

R 

EE:  ITS 

v yj ly  i m f\  rv  c. r.  l 

0 ERRORS  ASSEMBLY  COMPLETE 
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VORTEX  DASMR 


000®  HOURS 


I 


I 

i 


000002 

A 

1 B 

EOU 

2 

00O463 

A 

2 LHIJ 

EOU 

0462 

000005 

A 

3 TUNIT 

EOU 

c. 

0OOO06 

A 

4 DUO IT 

EOU 

6 

000004 

A 

5 AFC 3 

EOU 

4 

000000 

005111 

A 

6 START 

I AP 

OOOOOl 

054116 

A 

( 

STA 

FCB 

000002 

010006 

A 

2 

LDA 

DUNIT 

000003 

124007 

A 

9 

ADD 

REIJ+3 

000004 

054006 

A 

10 

STA 

REU+3 

O000O5 

010006 

A 

11 

LDA 

DUNIT 

00000b 

124020 

A 

12 

ADD 

READ+3 

000007 

054O17 

A 

13 

STA 

READ+3 

14  REU 

PEU 

FCB 

000010 

006505 

A 

00001 1 

000000 

E 

000013 

100000 

m 

000013 

OO140O 

A 

O0OO14 

000120 

R 

000015 

000000 

A 

O0PO16 

000000 

A 

000017 

03000  4 

A 

IS 

LDP 

AFCB 

OOOOcIO 

015003 

A 

16 

LDA 

3 , B 

000031 

0061 40 

A 

17 

SUE  I 

22 

0O0O33 

000036 

A 

O0OO23 

054077 

A 

13 

STA 

FCB+3 

15  READ 

READ 

FCB , , 

O00O34 

0O65O5 

A 

O0OO35 

000011 

E 

O00026 

1 00000 

A 

000037 

030000 

A 

000030 

000120 

R 

O0OO31 

oooooo 

A 

0OOO33 

000000 

A 

. 

000033 

002.000 

A 

20 

CALI. 

UR  I TRY 

000034 

000055 

R 

0 0 0 O 3 5 

001010 

A 

31 

JAZ 

EXIT 

000036 

000052 

R 

000037 

005010 

A 

32  F0R3 

LDA  1 

43 

000040 

000053 

A 

0OOO41 

1 40005 

A 

23 

SUP  ■ 

TUNIT 

000043 

050005 

A 

24 

STA 

TUNIT 

00004  3 

002003 

A 

35 

CALL 

UR  I TRY 

0000  4 4 

O0OO55 

R 

000O 45 

001010 

A 

36 

JAZ 

EXIT 

000O46 

P 00052’ 

R 

000047 

f 00 000 

A 

27 

HLT 

000060 

oo looo 

A 

3S 

IMP 

F0R3 

o floor,  i 

000037 

R 

29  EXIT 

CX1T 

0001V  .2 

001. '..0'. 

M 

000 O' . '■ 

oooooo 

1 

0000'. 4 

0002:00 

A 

27 


VORTEX  DASMP 


9009  HOURS 


000O55 

00000? 

A 

20  URITPV 

entr 

000056 

014006 

;A 

31 

LDA 

URITE+3 

000057 

150462 

A 

32 

ana 

LHU 

000060 

120005 

A 

3 3 

ADD 

TUNIT 

03OO61 

054033 

A 

34 

STA 

URITE+3 

35  URITF 

URITE 

RCB, 21 , ,4 

O0O062 

006305 

A 

0OOO63 

000025 

jr 

000064 

1 30000 

A 

000065 

0-10425 

A 

O0OO66 

0O012O 

R 

000067 

030000 

A 

000070 

000000 

A 

36 

STAT 

URITE, NEX" 

0O0O71 

006505 

A 

000072 

000000 

E 

000073 

000062 

R 

000074 

000 1 00 

R 

000075 

000130 

R 

OOO076 

030103 

R 

00O077 

O001OO 

R 

000100 

005001 

A 

37  NEXT 

TZA 

000101 

001000 

A 

38 

JMP* 

UP  I TRY 

000102 

103355 

R 

000103 

01 4306 

A 

39  EOT 

LDA 

SREC +3 

000104 

150462 

A 

40 

ANA 

LHU 

000105 

1 23005 

A 

41 

ADD 

TUNIT 

O0010G 

054OO3 

A 

42 

STA 

SREC+3 

43  SREC 

SREC 

FOB. 21 , , 1 

000107 

036505 

A 

0001 10 

000063 

E 

00011 1 

100300 

A 

0031 12 

012025 

A 

OO0113 

0OO12O 

R 

000114 

000000 

A 

000115 

0003O0 

A 

000116 

001000 

A 

44 

JMP* 

URITPV 

000117 

130055 

R 

45  PCE 

FOB 

2S60, BURR 

000120 

005000 

A 

000121 

0OO132 

R 

COO  122 

001400 

A 

060123 

0O3O30 

.A 

000124 

003000 

A 

000 1 £5 

000000 

A 

000 1 26 

000033 

A 

O00127 

030000 

A 

000130 

000000 

A 

0001 31 

009000 

A 

000  1 3,': 

46  BURR 

BSS 

2C60 

OOP0O0 

1< 

47 

END 

START 

F HTK’Y  1 

NAMf  5 

EXTERNAL  HAMMS 


VORTEX  DASHR 


000053  E VSEXEC  0001 10  E VSICC  000072  E V*I0ST 

SYMBOLS 

0OOOO4  A ARCB  00000a  A C 000138  R BURR 

000103  R EOT  000052  R EXIT  000150  P RCB 

000462  A LHU  000100  R NEXT  000024  R REhP 

O0O1O7  P SPEC  000000  P START  000005  A TUNIT 


000132  R BURR 
000120  P RCB 
000024  R RE f-I> 
000005  A TUNIT 


0000  HOURS 


000006  A DUNIT 
000037  P R0R3 
000010  R REU 
000053  E VSEXEC 


00O11O  E VSIOC  000072  E VtlOST  000062  R WRITE  0OO05S  P URITRV 
O ERRORS  ASSEMBLY  COMPLETE 
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VORTEX  DASMR 


0001  HOURS 


1 

EXT 

*SE, SEUf 

2 

NAtlE 

I*TPV 

000001 

A 

X 

EOU 

1 

000002 

A 

4 

S 

EOU 

2 

O00O0O 

01 4105 

A 

c; 

START 

LDA 

INU11 

000001 

024076 

A 

6 

LDB 

IA 

000002 

064110 

A 

7 

STB 

FCB+1 

000003 

0O1016 

A 

9 

JANZ 

NOSCEP 

00000  1 

000030 

P 

000005 

006020 

A 

9 

LDBI 

CALSEQ 

OOOOOt 

0OO1O3 

R 

10 

SCHED 

5,1,106. 

000007 

00650S 

A 

000010 

000000 

E 

0000 1 1 

01010S 

A 

000012 

143152 

A 

00O013 

142303 

A 

0000 1 4 

146720 

A 

0000 IS 

151323 

A 

0000 1 6 

014071 

A 

11 

LDA 

1ST 

O0OO17 

0O1004 

A 

12 

JAN 

DECRM 

000020 

0OOO66 

R 

13 

REU 

FCB , *0 

600021 

006505 

A 

000022 

000000 

E 

0O0O23 

1 00000 

A 

00002-1 

001412 

A 

00OO25 

000 1 1 £ 

P 

000026 

000000 

A 

000027 

000000 

A 

14 

NOSCED 

READ 

FCB , 10, 

000030 

006505 

A 

00OO31 

000022 

£ 

O0O032 

lOOOOO 

A 

000033 

030012 

A 

000034 

000112 

P 

O0OO3S 

oooooo 

A 

O0OO36 

000000 

A 

000037 

024040 

A 

15 

LDB 

IA 

000040 

016364 

A 

16 

LPA 

244,  E 

000041 

001010 

A 

17 

JAZ 

UIP 

000042 

000050 

P 

000043 

001002 

A 

18 

JAP 

DECRM- 1 

000044 

OOO06S 

R 

00004S 

044047 

A 

19 

I NR 

FCB  + 3 

©00046 

001 000 

A 

20 

Ji-tr 

NOSCED 

000047 

000030 

R 

000050 

0050.11 

A 

2 1 

UIP 

TBA 

O0OO61 

0061-10 

A 

22 

sum 

5 

000052 

000005 

A 

000063 

05404 l 

A 

23 

STA 

PCBi  3 

24 

URITE 

FOB, 10, 

000064 

006606 

A 

E 


VORTEX  DASMR 


0001  HOURS 


00O05S 

000031 

E 

000  OSS 

1OO00O 

A 

oeo ?s? 

030412 

A 

0000 so 

000112 

R 

000061 

OOOOO0 

A 

0000(52 

O0OOO0 

A 

00OO63 

001000 

A 

25 

JMP 

DECRM 

000064 

000066 

P 

00006S 

044027 

H 

26 

INR 

FCB+3 

000066* 

01 4020 

A 

2? 

DECRM 

LDA 

INUM 

O0O067 

005311 

A 

23 

I)AR 

000070 

054016 

A 

£9 

STtt 

INUM 

000071 

014016 

A 

30 

LDA 

1ST 

OO0072 

001000 

A 

31 

JMP 

000073 

000073 

000000 

A 

32 

1 4 TP  V 

BES 

0 

O0O074 

002000 

A 

33 

CALL 

SSE 

O00O7S 

000000 

E 

00OO76 

000002 

A 

34 

DATA 

2 

O00O77 

O00OO0 

A 

35 

IUNIN 

DATA 

0 

000100 

000000 

A 

36 

IA 

DATA 

0 

000101 

OO1O0O 

A 

37 

JMP 

START 

000102 

000000 

R 

C0O1O3 

10007? 

R 

3S 

CA'uSEO 

DATA 

( IUNII 

000104 

0001 10 

R 

0001  OS 

0001 1 1 

R 

COO  1 06 

OOOl 07 

R 

000107 

000000 

A 

39 

INUM 

DATA 

0 

000110 

000000 

A 

40 

1ST 

DATA 

0 

0001 1 1 

000012 

A 

41 

IUNOT 

DATA 

10 

42 

PCE 

F CB 

24S,  0 

0C0112 

000365 

A 

O0O11 3 

C 0OOOO 

A 

000114 

000000 

A 

0001  IS 

000000 

A 

0001 16 

000000 

H 

000117 

000000 

A 

000 1 20 

000000 

A 

000121 

000000 

A 

000122 

0OOOOO 

A 

O0O123 

oooooo 

A 

43 

END 

ENTRY  NAMES 
000073  R 1 4 T P V 
EXTERNAL.  NAMES 
000000  E 4 BUFF 
SYMBOLS 

00OOOO  F 4PURF 
000000  R DE  OPM 
OO0I07  R I MUM 
000030  P MC'SCFP 

000050  p mi r 

0 ERRORS  ASSE 


( IUNIN)*, 1ST, IUNOT, INUM 


000075  E SSE 

00007S  E *SE 
000  J 1 p.  R FCR 
0O011O  P 1ST 
000000  R START 
000001  A X 

mkl.y  foMRLtTr 


000010  F.  VSEXEC  0O005S  E VSIOC 


000002  A B 000103  P CALSEO 

000073  R IETPY  000100  R IA 

000077  R IUNIN  000111  R IUNOT 

000010  E VSEXEC  0000S5  E VSIOC 
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VORTEX  DASMR 


0000  HOURS 


1 

EXT 

S3E 

000001 

0 

2 

X 

EOU 

t 

000420 

A 

3 

ZERO 

EGU 

0420 

000423 

A 

4 

FOUR 

ECU 

0423 

000002 

A 

5 

B 

EGU 

2 

6 

EXT 

**BUFF 

0000  CO 

002000 

A 

? 

START 

CALL 

JfEUFF , < CALSEQ ) X , BUFR , D2560, ZERO , FOUR 

000001 

0OOOOO 

Z 

000002 

10003? 

R 

O0OO03 

0OOOS6 

R 

000004 

000055 

R 

O0O00S 

000420 

A 

O0OO06 

000423 

A 

00000? 

034030 

A 

8 

LDX 

CALSEQ+1 

000010 

O550O0 

A 

9 

ST0 

9,  X 

000011 

001004 

A 

10 

JAN 

ERROR 

000012 

000025' 

R 

000013 

014042 

A 

11 

LDA 

BUFR 

000014 

001010 

A 

12 

JAZ 

START 

0000 IS 

000000 

R 

0000 1 6 

064036 

0 

1 3 

STP 

NUDS 

14 

EXT 

DCMP 

000017 

002000 

0 

15 

CALL 

DC  MP , ( C0LEEO  + 2 ) * , (CALSEQ+3)*, BUFR. NUDS 

000020 

000000 

E 

O00021 

10004  1 

R 

000022 

100042 

R 

000023 

000056 

R 

000024 

000055 

R 

0 

16 

ERROR 

EXIT 

000026 

006505 

0 

* / 

0 O O 0 2 6 

OOOO0O 

Z 

00002V 

000200 

A 

. 

000020 

064002 

A 

17 

ENTRY 

STP 

t + 3 

00OO31 

001000 

A 

13 

JMP 

* + 3 

O00O32 

O0OO34 

R 

00OO33 

000045 

R 

19 

DATA 

TST 

000034 

002000 

0 

20 

CALL 

*SE 

00OO 36 

000000 

E 

000036 

000004 

A 

21 

DATA 

4 

000037 

000000 

A 

22 

CALSEQ 

DATA 

0, 0,0,0 

000040 

000000 

A 

000041 

oooooo 

A 

000042 

000000 

A 

000043 

00 1 000 

A 

£3 

JMP 

START 

000044 

000000 

ft 

000045 

R 

24 

TST 

EOU 

* 

25 

PUP 

4 

000045 

000051 

R 

26 

DnT  A 

t + 4 

000046 

000052 

ft 

26 

DATA 

T + 4 

00004? 

000063 

ft 

£6 

DATA 

* + 4 

000060 

00006.4 

P 

2i_. 

DATA 

* + 4 

000061 

00002V 

A 

27 

TST  1 

DATA 

23 .0,33,0 

000062 

000000 

A 

I 

VORTEX  DASMP  0000  HOURS 

O0OO53  00004 1.  0 
000054  00000 0 O 


O00O55 

005000 

A 

2S 

D2S60 

DATA 

2560 

0O0055 

R 

23 

NUDS 

EOU 

D2560 

- 

000056 

30 

BUFR 

BSS 

2560 

000030 

P 

31 

END 

ENTRY 

ENTRY  NAMES 
EXTERNAL  NAMES 

000001  E ft BUFF  000035  E «SE  000020  E DCMP  000026  E V1EXEC 
SVMBOLS 

OOOOOl  E f 3 BUFF  O0OO35  E *SE  000003  A B 00O056  R BUFR 

00O037  R CALSEC*  O00OSS  P 02560  O00O2O  E DCMP  000030  R ENTRV 

000O25  R ERROR  000423  A FOUR  O0O05S  R MUDS  000000  R START 

000045  R TST  000051  P TST 1 000026  E VtEXEC  000001  A X 

000420  A ZERO 

0 ERRORS  ASSEMBLY  COMPLETE 


I 

I 
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I 


VORT 

EX  DASMR  8BUFF  0000  HOURS 

1 

TITLE 

SBUFF 

2 

NAME 

% S BUFF 

3 

NhME 

LPECL 

4 

EXT 

*5E 

000421 

A 

r, 

ONE 

ECU 

0421 

000000 

000000 

A 

9 

8 f BUFF 

ENTR 

000001 

002000 

A 

*7 

CALI. 

SEE 

000002 

000009 

6 

000003 

000003 

A 

3 

D A T A 

5 

000004 

9 

UNIT 

B55 

1 

000005 

10 

WHERE 

BSS 

1 

C'0'0006 

11 

LENG 

BS5 

1 

000007 

12 

OP 

BSS 

1 

000010 

13 

MODE 

BSS 

1 

00001 1 

096017 

A 

14 

LDAE 

UHERE 

000012 

000006 

R 

, 

000013 

064047 

A 

15 

6TA 

LRFCL+1 

000014 

006017 

A 

16 

LDAE* 

LENG 

000016 

100006 

R 

00001 6 

064043 

A 

17 

5TA 

LPECL 

000017 

006017 

A 

18 

LDAE* 

MODE 

000020 

100010 

R 

000021 

004244 

H 

19 

LRLA 

4 

00OO22 

006127 

A 

20 

A DDE* 

OP 

000023 

1 00007 

R 

000O24 

OO4250 

A 

21 

LRLA 

3 

000026 

006127 

A 

22 

APDE* 

UNIT 

O0O026 

100004 

R 

O0OO27 

064003 

A 

23 

STA 

m 4 

000030 

R 

24 

WRITE 

ECU 

1 

26 

WRITE 

l.RECL,6,,l  WILL  BE  MODIFIED 

000030 

006505 

A 

000031 

000000 

E 

00O032 

100000 

A 

000033 

010405 

A 

00O034 

000062 

R 

000036 

000000 

A 

000036 

000000 

A 

000037 

006027 

A 

26 

LDBE 

WRITF+S 

000040 

'■>000  36 

P 

00004 1 

010421 

A 

27 

LDA 

ONE 

28 

5T  AT 

WRITE..  ERR.  FRR<  1 . ERP  + 2 , ERR  + 3 

000042 

006505 

A 

0 (i  0 0 4 y 

00 00 00 

E 

000044 

OO0030 

P 

00004  V. 

O00O6  i 

p 

000046 

000064 

R 

0OO047 

00006s 

R 

0000'  .,(> 

000066 

P 

0000'.  1 

00 1 000 

0 

29 

RE  TU* 

**BUf  F 

000062 

1 00000 

R 

000063 

006311 

A 

30 

ERR 

DAP 

000064 

0OS311 

A 

31 

DAR 

35 


I 


A 


VORTEX  BASMR  «BUFF  0000  HOURS 


000O55 

0O5311 

A 

32 

DAR 

000056 

O05311 

A 

33 

DAR 

00006? 

OO1O0O 

A 

34 

PETU* 

»*BUFF 

000060 

000061 

lOOOOO 

p 

35 

STA  BSS 

1 

36 

LRECL  FCE 

0,0,1 

000062 

000000 

A 

00OO63 

000000 

A 

000064 

000400 

A 

000065 

000000 

A 

00OO66 

000000 

A 

00006? 

000000 

A 

000070 

000000 

A 

000071 

000000 

A 

600072 

oooooo 

A 

000073 

000000 

A 

37 

END 

ENTRY  NAMES 

000OO0  P *IPUFF  000062  R LPECL 
EXTERNAL.  NAMES 

000002  E *SE  000031  E V*IOC  000043  E V*IOST 
SVMBOLS 

00O0O0  P *f BUFF  0O0002  E *SE  O000S3  R ERR  000006  R LENG 

000062  R LPECL  0O0O10  R MODE  000421  A ONE  00000?  R OP 

00OO61  R ST  A 000004  R UNIT  000031  E VilOC  0O0O43  E VtlOST 

000O05  R UHEPE  000030  R UPITE 

0 ERRORS  ASSEMBLY  COMPLETE 


VORTEX  DASMR  *BUFF 


0000  HOURS 


1 

NAME 

IUNPK 

000000 

006037 

A 

2 

START 

LDXii* 

NTOUD 

O0OO01 

100035 

•p 

O00OO2 

005017 

A 

3 

LDAE* 

IPK 

000303 

100034 

P 

000004 

124033 

A 

4 

ADD 

LLS 

O00O0S 

054010 

A 

5 

STA 

LLS1 

O00OO6 

014032 

A 

6 

LDA 

LASL 

000007 

005147 

A 

7 

SUBET 

IPK 

300010 

100034 

R 

000011 

054OO5 

A 

3 

STA 

LSL1 

000012 

006017 

A 

9 

LDAE* 

IA 

00O013 

1O0032 

R 

000014 

005344 

A 

10 

DXR 

DXR 

000015 

0O5OO2 

A 

1 1 

T2B 

000O16 

004540 

A 

12 

LLS1 

LLSR 

0 

000017 

004117 

A 

13 

LSL1 

ASRB 

IS 

000020 

006067 

A 

14 

STBEX 

IB 

000021 

100033 

R 

000022 

O44010 

A 

15 

I NR 

IB 

000023 

001046 

A 

16 

JXNZ 

DXR 

000024 

000014 

R 

O00O25 

001000 

A 

17 

JMP 

* 

0OOO26 

000025 

P 

000026 

IB 

IUNPK 

EES 

0 

19 

EXT 

SSE 

O0OO27 

0O20OO 

A 

20 

CALL 

*SE 

000030 

00.3000 

E 

000031 

000004 

A 

21 

DATA 

4 

000032 

000000 

A 

22 

IA 

DATA 

0 

00OO33 

0OO000 

A 

23 

IB 

DATA 

0 

00O034 

000000 

A 

24 

IPK 

DATA 

0 

000035 

000000 

A 

25 

NTOUD 

DATA 

0 

000036 

001000 

A 

26 

JMP 

START 

000037 

0OOO0O 

R 

000040 

004540 

A 

27 

LLS 

LLSR 

0 

000041 

0041 17 

A 

2S 

LASL 

ASPB 

15 

29 

END 

ENTPY  NAMES 
000026  P IUNPK 
EXTERNAL  NAMES ' 

O0OO3O  E *SE 
SYMBOLS 

O00O30  E *SE  000014  R DXR  000032  R IA  000033  R IB 

O0O034  R 1 F’K  000026  R IUNPK.  00004  1 R LASL  000040  R LLS 

O0001S  P 1.1. A 1 00OO17  R LSI.  1.  O0PO3G  P NTOUD  000000  R START 

0 ERRORS  ASSEMBLY  COMPLETE 
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VORTEX  DASMR  *BU FF 


0000  HOURS 


I' 


f 


i 

J 

> 


i,  i 


1 1 


< 


L 


i 


2 

NAME 

I SHF 

000001 

A 

3 X 

EOU 

1 

000000 

006037 

A 

4 START 

LDXE* 

N 

000001 

100022 

R 

000002 

00S344 

A 

S DXR 

DXP 

000003 

COSE  1 S 

A 

6 

LDAE* 

I A , X , 0200 

003004 

1 00020 

P 

000O0S 

004244 

A 

7 

LRLA 

4 

000006 

004304 

A 

s 

ASPA 

4 

000007 

0062S5 

A 

9 

STAE* 

IB, X, 0200 

0000 1 0 

100021 

R 

000011 

001046 

A 

10 

JXNZ 

DXR 

000012 

000302 

P 

000013 

001030 

A 

11 

JMP 

* 

000014 

000014 

000013 

R 

12  I SHF 

BES 

0 

13 

EXT 

*SE 

O0O01S 

002000 

A 

14 

CALL 

SSE 

000016 

OOOO0O 

E 

000017 

000003 

A 

IS 

DATA 

3 

000020 

16  I A 

BSS 

1 

000021 

17  IB 

BSS 

1 

000022 

IS  N 

BSS 

1 

000023 

O01OOO 

A 

19 

•JMP 

START 

000024 

O0O0OO 

R 

20 

END 

ENTRY  NAMES 
000014  P I SHE 
EXTERNAL  NAMES 
000016  E tSE 
SYMBOLS 

000016  E *SE  000002  R DXR  0O0020  R IA  000021  R IB 

000014  R I SHF  003022  P N 000000  P START  000001  A X 

0 ERRORS  ASSEMBLY  COMPLETE 
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VORTEX  DASMR  *BUFF  0000  HOURS 


1 

2 

NAME 

S*BF1 

3 

EXT 

SSE 

000000 

014044 

0 

A 

START 

LBA 

FRST 

O0O001 

001016 

0 

5 

JAN  7. 

URITE 

000*502 

000025 

R 

00O0O3 

044041 

A 

6 

INR 

FRST 

000004 

00602? 

A 

7 

LDBE* 

UNIT 

00  00. 35 

100041 

R 

O00006 

0O5021 

A 

8 

TBA 

000007 

124011 

A 

9 

ADD 

REU+3 

000010 

054010 

A 

10 

STA 

REU+3 

000011 

OO5021 

A 

11 

TBA 

000012 

124015 

A 

12 

ADD 

WRITE+3 

000013 

054014 

A 

13 

STA 

URITE+3 

000014 

0 1 4035 

A 

14 

LDA 

IB 

000015 

054031 

A 

15 

STA 

FCB  + 1 

16 

RED 

REU 

FCB 

000016 

006505 

A 

000017 

000000 

E 

000020 

100000 

A 

000021 

001400 

A 

000022 

000046 

R 

000023 

000000 

A 

O00O24 

000000 

A 

1? 

WRITE 

URITE 

FCB, , ,3 

00OO25 

006505 

A 

000026 

O0OO17 

E 

00002? 

100000 

A 

00OO3O 

030400 

A 

O00O31 

000046 

P 

0O0032 

000000 

A 

000O33 

000000 

A 

000034 

0O100O 

A 

12 

JMP 

* 

000035 

000034 

P 

O00O35 

19 

*tBFl 

BES 

O 

000036 

002000 

A 

20 

CALL 

*SE 

000037 

0O0O0O 

E 

000040 

O0OOO2 

A 

21 

DATA 

2 

000041 

000000 

A 

22 

UNIT 

DATA 

0 

000042 

000000 

A 

23 

IB 

DATA 

0 

000043 

001000 

A 

24 

JMP 

START 

000044 

000000 

P 

000045 

000000 

A 

25 

FRST 

DATA 

0 

26 

FCB 

FCB 

245,0,1 

000040 

000365 

A 

000047 

0O00O0 

A 

000050 

000400 

A 

000051 

000000 

A 

000052 

000000 

A 

000053 

000000 

A 

000054 

000000 

A 

O00O55 

000000 

A 

VORTEX  DASMR  SBUFF  0000  HOURS 

000056  0O0O00  A 

OOOOS?  O00O0O  a 

2?  END 

ENTRY  NAMES 
000035  R t-  * E p 1 
EXTERNAL  NAMES 

O0OO37  E *SE  00002S  E V*IOC 
SYMBOLS 

000035  R **BF1  000037  E *SE  000046  R FCB  000045  R FRST 

000042  R IB  000016  P PEU  O00000  R START  000041  R UNIT 

000026  E V*IOC  00OO25  P WRITE 

0 ERRORS  ASSEMBLY  COMPLETE 


0000  HOURS 


VORTEX  FTN  TV 


3 SUBROUTINE  I BCMP ( I UN  IT , INUrt , NTMP , I A , NN ) 

4 DIMENSION  IAC  1 ) , IBf.245) 

5 IF ( IAC 1 > , NE . 1 )G0  TO  39 

6 I NUN* INUM+1 

7 CALL  ISHFC IAC 4) , IAC4)  . 1 ) 

8 1ST  »IA<4)+1 

9 NUBS-240 

10  IFCIST.LT.0iGO  TO  1 

11  3 IFCNTMF  . NE  . GOCALL  *» BP  1 C I UN  IT , I B > 

12  NB  *240 

13  I OUT  * 1 

14  NTMP-0 

15  I F ( 1ST. GT. 0 ) GO  TO  4 

16  CALL  ISHFC IACg) , IB.NUDS) 

17  GO  TO  5 

18  4 Mil*  16 

19  DO  11  1*1 .3 

20  LL*S0*I-49 

21  11  CALL  JDCMPC IAC 1+4 ) , I ACMM)  , IBCLL) . S0, MM) 

22  DO  12  1*1,3 

£3  LL  *1041  + 141 

24  12  CALL  JDCMPCIACI+7)  , IA.'MM)  , IBCLL)  , 10, MM) 

25  Call  ISHFC IACMM) . IBC131 > . 9) 

26  MM-MM+9 

£7  CALL  JDCMPC IAC 11 ), I ACMM) , IBC 190) , £3, MM) 

28  CALL  JDCMPC IAC 12) , IACMM) , IBC213) . 1S.MM) 

29  CALL  ISHFC IA(MM) . IBC223) , 13) 

30  NBEG-13 

31  6 IBCMB+l ) -IAC4 ) 

32  IBC NB+2) * IACNBEG) 

33  IBCNB+3) -IACNBFG+1 ) 

34  CALL  ISHFC IAi NBFG+2) , IBCNB+4) , 1 ) 

35  I8CNB+5) »NUDS 

36  IF CNTMP . NE . 0 > IBCNB) »-l 

37  IF CNTMP , EO . 48) I0UT*1 

38  IF ( I OUT . EO . 0 ) RETURN 

39  98  CALL**EF1( I UNIT, IB) 

40  NTMP-0 

41  RETURN 

42  99  IF  C HTMP . EO . 0 ) RETURN 

43  IF  < J A C NM+ 1 ) . EO . 409S ) GO  TO  98 

44  RETURN 

45  1 NUID5-  IAC  3) -7 

46  If  OIMUC.Hl  . com  TO  3 

47  I OUT  * O 

4f:  NIUC.  5 

49  NTMP'-NTMC  + l 

50  NP.*f:45-5TNTMP  • 

51  GO  TO  6 

52  END 
ENTPV7C0MM0N  BLOCK.  NnMFS 
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VORTEX  FTN  IV  0000  HOURS 


00144? 

R 

IDCMP 

EXTERNAL 

NAMES 

000002 

E 

*se 

001131 

E 

I SHE 

001235 

£ 

SSEF1 

001002 

E 

JDCMP 

000733 

C 

*L0 

SYMBOL  TABLE 

0O136? 

R 

000001 

001444 

R 

000002 

1O0OO4 

R 

I UN  I T 

1O00O5 

R 

I HUM 

100006 

P 

NTMP 

100007 

P 

I A 

100010 

R 

NN 

0OOO02 

r 

f SE 

O0O013 

R 

IB 

001446 

p 

000365 

001237 

R 

99 

0O1131 

E 

I SHF 

001366 

R 

000003 

O01366 

R 

*1 

091356 

R 

0* 

001370 

R 

1ST 

001372 

P 

NUDS 

O01371 

R 

000360 

0O337 3 

P 

000000 

001 306 

R 

1 

000467 

R 

3 

001225 

E 

S*EF  1 

001374 

R 

NB 

001375 

R 

I OUT 

000S52 

R 

4 

O01376 

R 

000007 

001033 

R 

6 

001400 

P 

MM 

00] 377 

F! 

000020 

000670 

R 

I 1 

00 1401 

F' 

I 

00 1 404 

R 

l.L 

001402 

P 

000062 

O014O3 

R 

000061 

001002 

E 

JDCMP 

0O1405 

P 

173777 

031 406 

P 

f 1 0 

00] 407 

P 

OOOO 1 2 

0O 14 10 

P 

>1  1 

0O07O  3 

f; 

*D0 

000646 

F' 

1 2 

0014)1 

F> 

OO0O  i 

00 14  11! 

P 

000,3  1 • , 

00 1113 

P 

000001 . 

001 415 

p 

0000 1 1 

42 


0000  HOURS 


VORTEX  PTN  IV 


031414  R 000277 
09 1420  R 000013 
091417  R 000027 
001416  P 000310 
00142a  p 009017 
001421  P 000337 
001424  R 0000  IS 
001423  R O0035S 
001425  R NF.EG 
001426  R 900013 
001427  P OOO014 
001430  R *12 
901431  P *13 
001432  E *1  4 
001433  R 0000 1 S 
001434  R *1  S 
001435  fi  *16 
001436  R *17 
00 14 3?  P 0O0M16 
001440  P *13 
901441  P 000017 
001442  R 000069 
O01224  R 38 

00144?  p 007777 
001445  P OOO0O5 
0 ERRORS  COMPILATION  COMPLETE 


VORTEX  FTN  IV 


0000  hours 


1 

SUBROUTINE  JDCMRi IPX, IP, IB, N, MM ) 

2 

DIMENSION  IB<N),IAf2) 

3 

IF (IPX ,N£. 12100  70  1 

4 

CALL  ISHFt  IA,  IB, II' 

5 

MM-MM+N 

6 

RETURN 

7 

1 

IF  t IRK . NE . 0 >60  TO  2 

s 

CALI.  ISHFt  In,  IB.  1 > 

'3 

DO  22  I » 2 , N 

10 

22 

IBCn»IB(  1 ) 

11 

MM-MM+1 

12 

RETURN 

13 

2 

CALL  ISHFt Ifi, IB, 1 ) 

14 

NTOUD* 127IPK 

15 

L-a 

16 

NN*N-1 

17 

DO  3 1*1, NN , NTOUD 

1? 

CALL  IUNPK <IA(L),IB(I*1>, I PK, NTOUD) 

19 

L*L  + 1 

20 

3 

CONTINUE 

21 

MM  * MM+L 

22 

DO  4 I * 2 , N 

23 

4 

IP(I)*IE(I)+IB(I-1) 

24 

PETURN 

SC 

END 

ENTP 73 COMMON 

BLOCK  NAMES 

00O341  R 

JDCMP 

EXTERNAL  NH»iES 
000002  r *SE 
0001 -IS  b'  I SHF 
0003 IS  E *D0 
00028*.  E IUNPK 
SVMEOl.  TOliLt: 

000326  P 000001 
OO0327  R 00O002 
100004  R IRK 
100OOS  P 10 

100006  R IB 

100007  R N 

1O0010  R MM 

000002  F.  *SE 
000324  R 000014 
000311  R 0S 

OO0OE4  P 1 

OOOl 4S  E I SHF 
00O32S  R 000000 
O00144  R S 

0OO1O1  R 22 

000330  P I 

OO0331  R 177777 
000332  P *1 
0003  IS  E SDO 
000333  R NTOUD 

VORTEX.  FTN  IV  0000  HOURS 

000334  R L 

00033S  P NN 

000237-  R 3 

000225  F.  IUNPK 
0OO736  P *1  0 
000264  P 4 

000237  P 17777 6 
000V 40  R tl  t 
0 ERRORS  COMPILATION  COMPLETE 
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Each  physical  record  consists  of  2560  or  fewer  12-bit  words,  (512  CDC  6600 
words),  and  contains  as  many  logical  records  as  can  be  completely  contained 
in  2559  words.  All  logical  records  contain  3 control  words,  except  for 
type  -1. 

Word  1 * type  of  record 
Word  2 = subtype 

Word  3 * number  of  12-bit  words  this  logical  record 
(including  control  words). 

At  present,  three  types  of  records  have  been  reserved. 

Type  0 * header  record,  not  yet  defined 
Type  1 = frame  data 

Type  -1  = end  of  physical  record  (12-bit,  2's  complement  -1). 

For  type  1 records,  subtype  0 has  been  defined. 

Word  1 = 1 
Word  2=0 

Word  3 = N = number  of  words 

1 - good  frame 

Word  k = Status  = -1  - error  frame 
-2  - timeout 

The  rest  of  the  frame  is  as  follows  for  the  various  frame  statuses. 

Status 

-2  -1  -1 
"received"  time  seconds 

" " milli seconds 

Buffer  contents 

partial  frame  240  data  words  received  in 

received,  if  arranged  form  by  instrument 

any.  as  given  in  AFCRL-TR-75_0588 . 


Word  5 
Word  6 
Word,  7 
Word  8-N 


Data  and  buffer  contents  are  two's  complement 


11  10 


1 0 


s 


data 


e 


Status  1 

Word  5 “ packing  factor  SCX 


Word  6 - " " SCY 

Word  7 - " " SC2 

Word  8 - " " FFX 

Word  9 - " " FFY 

Word  10  - " " FFZ 
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Word  11  - packing  factor  Analog  spares 
Word  12  - " " Digital  spares 

The  rest  of  the  frame  is  then  the  same  as  words  5-N  of  status  -1 
frames  except  that  relevant  data  is  compressed  according  to  the 
scheme  detailed  below. 

Packing  Scheme 

If  a particular  instrument  gives  in  readings,  X^, , X^, 

m-1  AX.  are  calculated  where  AX.  • Xj+1  -X..  If  all  the  AX. 
are  0,  the  packing  factor  K is  chosen  to  be  0.  Otherwise  K is 
the  minimum  of  2,  3,  **,  6,  12  such  that  -2*  ' < AXj  < 2K  ' -1  for 
all  i * 1 , . . . , m= 1 . 

When  K-12,  actual  data  is  stored.  When  Kj<12,  the  packed  instrument 
format  is  as  follows,  where  l * 12/K 


^ for  K M 


It  should  be  noted  that  X,  and  X are  stored  with  error  bits  (which  are 

1 m 

zero  because  we  only  pack  error-free  frames)  but  that  the  AX.  are  data 
deltas  only,  in  two's  complement  form. 

Therefore,  after  the  AX.'s  are  expanded  to  full  words,  X.  = Xj_^+2AX._^ 
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2.  MATHEMATICAL  ANALYSIS 

The  main  activity  of  the  past  year  has  been  writing  programs 
to  handle  large  amounts  of  data. 


(a)  MULTIPLE  CORRELATION 

A multiple  correlation  program,  MLTCOR,  (Program  1)  was 
written  that  yields  the  partial  regression  coefficients 
between  ground  measured  geophysical  parameters  such  as 
Kp,  Ap  and  AE  and  various  solar  wind  parameters  such  as 
bulk  velocity,  magnetic  field,  and  the  direction  of  the 
magnetic  field.  Using  the  derived  coefficients  predictions 
were  made  of  Kp,  Ap,  and  AE  and  these  could  be  plotted 
and  compared  with  the  observed  values. 
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Program  1 


GELKCt  CK1 10 000 » T150  . 

-ATTACH,T4FE1,C0R0TX3213,ID=BELLEW,MR=1. 
ATTACH, TAP;2,IMFX3218, IO=BELL£W,MR=i • 

ATTACH, PEN, ONLINEPEN, MR=1 . 

FTN,SL,PL=77777. 

•LOSET  (LIB=PEN) - 

LGO. 


2261 


BEL  LE  H 


10 

11 

12 

•13 

14 

15 
. 16 


l 

t j 


PROGRAM  ML TCOR < INPUT, OUTPUT, T APE  1, TAPE  2, TAPE  3, TAPE 6=OUTPU 
♦ T) 

CATALOGED  AS  MLTC0RX3213,CY=5 

DIMENSION  0 (6  > , SB  (6  ) • 

DIMENSION  PROGIO ( 3) 

COMKON/SU9/ 3 , SB , SY , BNOT 

COHMON/OAT/LAG, JMPOT 

OAT A PR0GI0/6RBELLEW, 4R2261,2RB9/ 

CALL  PLTIC3  (^RGGIO, 23 C . 0 til .3, 0 .75) 

F=2,5  • 

TOL=.00Ql 

N=6  - - 

H=1  00 
SIGNs-1 
JKPDT  = 0 

DO  2 J = l,  2 - ~~  - - 

LAG=1 

DO  1 1 = 1,1  

CALL  SETUP(N,M,STSN) 

- CALL  SU91  (F,TOL,M)  - - - 

LAG=LAG+i 

- JMPOT  = JMPDT+l 

CALL  ENOPLT 

--STOP-  — 

END 


13 

13 

11 

12 

13 

14 

15 
15 

17 

18 
19 
22 
21 
22- 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 
- 34 

35 

36 
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SUBROUTINE  SETUP (NtMtSIGN)  3 7 

• 33 

SETUP  REA03  0 AT  A FROM  TAPEl  ANO  CONSTRUCTS  39 

A MATRIX  OF  THE  CROSS  CORRELATION  COEFFICIENTS  43 

41 

OIMENSION  A (&,6) ,Y( 600 ,6) , YB(o> ,SIG(5> ,Rf 6,6) - 42 

CGMMON/SET/Y3,SIG,P,NM,OF,Y  43 

COMMON/D AT/L A G»  JMPD  T - 44 

CCMMON/MUL/KNT  45 

REWIND  3 4c 

PRINT  ICC  47 

FOR  MAT  ( 1H1 ) • - — 48 

FORMATO  49 

- PRINT  3 00  , SIGN, LAG  - - 50 

FORMAT (/*  SIGN=*,F5.1,*  IAG=*,I2///)  51 

N NUMBER  OF  DEPENDENT  AND  INDEPENDENT  VARIABLES  52 

M NUMBER  OF  OBSERVATIONS  . 53 

- HSAV=M  - — - - - . 54 

CALL  OATINIM)  55 

- PI=2.0*ACOS(0.G> r — 56 

FP=1.0/(4*PI)  5 7 

REWIND  3 ■-  - - - 58 

NM= N-l  59 

• - OF=M-l - - — - " ' . 63 

KNT=0  61 

— 00  1 J=1,M  _ 62 

READ( 3)  IDATE,8,BX,3Y,8Z,SG8X,SG3Y,SGBZ»  63 

«-VTH,RHO,V,VX,  VY,VZ,  AE  — 64 

PRINT  75,  8 , 8X  ,BY  , BZ  , VTH » RHO  tV,VX»VY,VZ»AE  65 

FORMAT  (4F1C  .4/7F10. 4)  - • • 66 

IF(SIGN*BZ.LT.C.O>  GO  TO  1 67 

- KNT=KNT  + 1 .-63 

YCKNT  ,l)  = RHO*V*V  69 

— -SX*FP*((BY* 3Y*3Z»8Z>*VX-(VY*3Y+VZ*3Z)*9X)  - 70 

SY=FP*((3X*3X+BZ*BZ)*VY«(VX*3X+VZ*3Z>*BY)  ' -71 

- SZ=FPM  (3X*3X*6Y*BY)*VZ-(VX*3X  + VY*3Y)*3Z>  72 

S=SQR.T(SX*SX+SY*SY*SZ*SZ)  73 

--PRINT- 11, S,SX,SY,SZ  - - 74 

YCKNT,  2)=S  . 75 

Y CKNT , 2 ) = BZ  - — 76 

YCKNT,  3>  = V*3Z  77 

Y CKNT , 5 ) = SORT (SGBX* SGBX+SG3Y*SG3Y+SG3Z*3S 3Z)  - 78 

Y CKNT ,4)=BZ*8Z*V  79 

Y CKNT, 2) = YCKNT, 3)  *Y CKNT, 5 1 - 81 

YCKNT, 6)  = AE  81 

CONTINUE  82 

M=KNT  83 

- PRINT  15,  N,  M - 84 

FORMAT  C/2110/)  85 

PRINT  11,  C (Y(J,I) ,1=1, N> , J=1,H>  , 85 

PRINT  200  87 
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SUBROUTINE  SUei<F,TOL,M)  129 

•••  * 133 

DIMENSION  0 Co ) f S3 (6 ) 131 

DIMENSION  Y3C6)  , SIG  (6)  ,R(6,6> : 132 

DIMENSION  Y C6CC  ,6)  133 

- COMHON/SUB/B,Se,SY,BNOT --  — - — 134 

COMMON/SET/YB,SIG,P.,NM,OF,Y  135 

- • ........  136 

PRINT  200  137 

FORMAT  (>  - - - 138 

CLRL=Q  139 

- F 1= F • - — - 14  j 

N=NM*1  141 

-F2=F 142  - 

PRINT  940,F1,F2  ls3 

- FORMAT  {*  F1,F2  =*1P2E12.3) - 144 

VMIN=1O**10  145 

VMA  X=0 . 0 1 146 

NMAX=  0 1h7 

~NMIN=0  - ■—  • — - - 143 

SY=SIG(N)  *SQRT(R(N,N)/OF)  149 

PRINT  120, SY — - - 150 

FORMATC*  SY=*F10.5)  151 

DO  10  J=1 , NM  • - *- 152 

B C J 1 = 0 153 

r 154 

1=1  . 155 

-CONTINUE 156  - 

PRINT*, R(I, I) ,1,1  157 

— * — • 159  • 

THIS  TEST  REDUCES  THE  POSSIBILITY  OF  DEGENERACY  159 

WHEN  AN  INDEPENDENT  VARIABLE  IS  APPROXIMATELY  A - - •■-•160 

LINEAR  COMBINATION  OF  OTHER  INDEPENDENT  VARIABLES.  161 

- IFCRCI,I)  .LE.TOL)  GO  TO  3 - - -.162 

* 163 

LINEAR  COMBINATION  OF  OTHER  INDEPENDENT  VARIABLES  164 

VI=R(I,N) *R(N,I)/R(I,I)  165 

PRINT*  ,VI,  VMAX,I  -- — - - 166 

IF(VI.EQ.O)  GO  TO  3 -167 

■ IF(VI.LT.0. 01  GO  TO  1 - 163  •• 

IF( VI . LE.  VMAX)  GO  TO  3 169 

VMA X* VI  - - 170 

NMAX=I  171 

- - PRINT*,  VMAX, NMAX  - - - - - 172 

GO  TO  3 173 

- B(I1=SIG(N)/SIG(I)*R(I,N)  17h 

PRINT*, 8(1), I 175 

- SB(I>=SY/SIG<I)*SQRT<1.0/R(I,I>)  — 176  - 

SB(I)=SY/SIG(I)*SQRT(R(I,I) ) 177 

PRINT  430,9(1) ,SB(I), I - * - - 178 

FORMAT (1P2£12«3,0PI10)  179 


I 


IF! AQSCVI ) .GE.A9S (VMIN) ) GO  TO  3 
\fHIH=VI 
,nhih=i 

IF(I.iQ.NH)  GO  TO  100  

1 = 1*1 

GO  TO  4 - 

SUM=0 

00  101  J=1,I 
SUM=SUM*3 (J)*Y9( J) 

3N0T=Y3(N) -SUM 

IF(CLRL,NE. v)  CALL  COMLC(B) 

FIN=A9S(VHIM)*0F/R(N,N) 

PPINT  710 ,FIN,F2, BNOT 

FORMAT (*  FINt  F2 , BNOT  =*1P3E12.3) 

IF(FIN.GE.F2)  GO  TO  6 

K=NMIN  - - - -- 

PRINT  620  , K 

FORMAK*  - K=NMIN=*I 2 ) • - ' ■ 

DF=OF«-l 

GO  TO  5 - - 

CONTINUE 

FOUTMMAX*  (OF-l.O)  / CR  ( N , N)  - \/MAX) 

PRINT  7 11,F0UT,F1 

FORMAT  FOUT.Fl  =*1P2E12.3> 

CONTINUE 

- IF ( MAX* (DF-1.0)/(R(N, N) -VMAX ) .LE.F1)  RETURN  - 
K=NMAX 

PRINT  123, K 

FORMAT (*  K=NMAX=*I2) 

- OF=OF-l  • --  - - - 

CALL  NUMA T { K, N) 

CLRL=1  -•  - 

GO  TO  7 


END 


18  j 
181 
182 
163 

184 

185 
185 
167 
189 
189 
193 

191 

192 

• 193 
19h 
195 
195 

197 

198 

199 
20u 
201 
202 
203 
2C4 

- 205 

205 

207 

' 20  3 

♦ 209 
210 
211 
212 

213 

214 
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SUBROUTINE  NUMAT ( KrN)  215 

OIMENSION  Y8<6>  ,SIG(6)  ,R(6,6),0!6,6I  ' • 216 

OIMENSION  Y(6CC,6)  217 

C0MM0NI/SET/Y3,SIG,R,NM,DF,Y 218 

PPINT  20C  219 

PRINT  50,  <(R(I,J>  ,J=1,N)  ,I  = i,N>  - — - — 223 

F0»MATC1P6E12.3)  221 

FORMAT!)  222 

00  1 1=1, N 223 

00  1 J=l, N 22* 

IF(I.EC.K.ANO.J.EQ.K)  GO  TO  10  225 

IF(I.cQ.K.ANO.J.NE.K)  GO  TO  9 226 

IFCI.NE.K.AND.J.EQ.  <)  GO  TO  8 227 

0(I,J)  = (P(I,J)*R(K,tO-R(I,K)*R(K,J)>/R<K,K>  - — ■ -228 

GO  TO  1 229 

0(I,J)=-R(I,K)/R(K,K)  - - 230 

GO  TO  1 * 231 

• 0(I,J)=R(  K,  J)/P(K,K) - - 232 

GO  TO  1 233 

0(I,J)  = 1.0/R(I,J1  5 234 

CONTINUE  235 

00  2 1 = 1,  N - - - • . • 236 

00  2 J=  1 , N « 237 

R(I,J)  = 0(I, J)  --  - - • 238 

PRINT  200  239 

-PRINT  50, (!R(I,J3 ,J  = 1,N) , I = i,N)  240 

RETURN  241 

END - ^ -242 
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SU9 ROUTINE  OATIN(M)  243 

OIHENSION  X tl-)  , IX(  1-.)  , 10(14)  * 244 

COMMON/OAT/LAG,  JMOOT  2w5 

IF< JMPOT.NE.QI  GO  TO  1C  ' ■ 246 

REWINO  1 247 

REWINO  2 . - 243 

REWINO  3 2-9 

LAGL=LAG  25 j 

LAGM=1  251 

LLAG=L AG- 1 - 25? 

KCUNT=0  253 

RSAD(l)  IOATElt (X(I),I=i,13)  254 

IFCEOF(l))  5,2  255 

REAO( 2 1 I0ATE2,IX11  256 

IF (EOF (2) ) 7,9  25 7 

XC14)=IX11  ■ " 253 

IFU0ATE1.GT.IDATE2)  GO  TO  2 . 259 

IFU0ATE2.EQ.IDATE1)  GO  TO  4 - - 263 

REAO(l)  I 0ATE1,(X(I), 1=1,13)  * 261 

IFCIDATE1 .GT.IDATE2)  GO  TO  2 r - 262 

IF< EOF ( II  ) 5,3  263 

j CONTINUE  •--•••  264 

IF(LAG.EQ.O)  GO  TO  12  . ..  265 

IFCLAGL.EQ.LAG.CR.LAG.EQ.l)  GO  TO  18  266 

j 00  J 7 1 = 1 , LLAG  267 

-ID(  I)  =ID(  1 + 1)  - -263 

IXtl) =IX  < I ♦ 1)  . 269 

CONTINUE - - -273 

I=LAG  271 

-GO  TO  22  - - - - 272 

1=1  273 

REAO(  2 ) I D d ) , IX  ( I)  - — - - - 2?t 

1X1 1=1 X ( I ) 275 

••  1 00 1=10 AT tl-(  10 AT £1/100*100 — - ---.276 

ID02=ID(I)-(IO(I)/100)*10C  277 

IDI  F=  I002-I  00 1 - 278 

IFdOIF.EQ.LAG.  OR.  IOIF.EQ.  LAG-24)  GO  TO  19  279 

• IFd.EQ.LAG.  ANP.LAGL.EQ.LAG)  GO  TO  19  - 280 

IF(IOd)-IDATEl.EQ.LAG)  GO  TO  19  231 

GO  TO  22  283 

L AGH=L  AG  — * — — ■ — .....  . 2 84 

X(14)=IX11  285 

WRI  TE  ( 3 ) I0ATE1 , <X(  I)  , 1=1*14) 286 

K0UNT=K0UNT4i  287 

PRINT  54,  IDATE1,  CXtl)  ,1=1,7) 288 

FORMAT(I10,7F10.3)  289 

-PRINT  56,  <Xd),  1 = 6,  14) -290 

FORMAT  (7F1C  .3)  291 

- IFTKOUNT. EQ.M)  RETURN  --  *292 

LAG  1=1  293 


IFflAG.EQ.O)  GO  TO  1 294 

- 1*1  - 295 

REAOfl)  I DA  Til, (X(J),J=1,13)  295 

• I F(  EOF  ( 1 ) ) 5,55  - - •* • 297- 

CONTINUE  293 

IFdDATEl.EQ.IDd>>  GO  TO  4 299 

1=1+1  3C3 

IF(I.LE.LAG>  GO  TO  55  . ..  3il 

LAGM=1  30? 

L A5L=l  AG  - — — — - — • ■ ■ • - 30  3 

GO  TO  2 30, 

PRINT  6 - 305 

FORMAT <*  EOF (1)  *>  306 

- RETURN  - 307 

PRINT  8 308 

FORMAT  (*  - EOF  (2)  *1 - 309 

RETURN  3io 

ENO  *•  * — _ — — ....  — • — — . — . ~ — . — — 311 

SUBROUTINE  COMLC(B)  312 

313 

THIS  SU°ROUTINE  CALCULATES  THE  TOTAL  CORRELATION  31<* 

- COEFFICIENT  OF  THE  OBSERVED  DEPENDENT  VARIA8LE  V,  315 

AND  THE  PRWDICTED  DEPENDENT  VARIABLE  YJ.  316 

317 

DIMENSION  Y8<6> ,SIG (6> ,R<6,6> ,Y (60? ,6>  318 

DIMENSION  B (6 ) 319 

DIMENSION  PY(12C» ,PX(103>  323 

C0MM0N/SET/Y3,SIG,R,NM,0F,Y  321 

CCMMON/MUL/M  322 

NrNM+i  - - - 323 

CYJY=0  324 

SIGYJsO  " - - 325 

OO  2 1=1, M 326 

YJ=0  327 

DO  1 J=1,NM  . 328 

YJ=YJ+8(J>*Y(I,  J)  329 

PY(  I)  =Y  (I , N)  330 

PX(I»=YJ  - — - - - 331 

CYJY=CYJY+YJ*Y(I,N)  332 

SIGYJ=SIG  YJ+YJ+YJ  333 

SIGYJ=SQ?.T  (SIGYJ)  33-, 

RML=CYJY/(SIGYJ*SIG(N>  > — - * 335 

P=0  336 

DO  3 1 = 1, NM 337 

IF(Bd).NE.O)  P=f>  + 1 338 

- CONTINUE  - - - 339 

RR=RML*RML  340 

FTEST=RR/ (1-RR) MM-P-1)/P  - 341 

PRINT*, M RMUL, FTEST,M, P=M , RML, FTEST, M,P  342 

CALL  SPLOT(PX,FY,M, 99 , XM AX, YMA X , Xu , Y0>  343 

CALL  PLOTTER(PX, PY, -1000. 0,10GC. 0,5. 0,0. 0,6, 11,1, Ml  344 

RETURN  345 

END  '*46 
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SU3P0UTINE  SPLOTIX, Y,N,ISCAL,XHAX,YMAX,XQ  .YC) 

0 1 HENS  I ON  X(1),Y(1)  ,G<101,61> 

OATA  BLANK/ 1H  /, AXIS/1H./,»0INT/1H*/ 

• ** 

* X(I)  = ARRAY  FOR  X-AXIS 

* Y ( I ) = ARRAY  FOR  Y-AXIS  • — • - - — ■ - 

* Ns  NUMPER  OF  POINTS 

* X0,Y0=  ORIGIN  OF  Gr APH 

* OX,OY=  SCALE  FACTCC  MULTIPLYING  POINT  VALUES 

* **NOTE**TOO  HIGH  OF  TOO  LOW  X,Y  WILL  CAUSE  ERROR  IF  ISCAL.NE. 

♦ 99 

* IF  ISCAL=99  THEN  DO  AUTO-SCALING 

IFCISCAL.NE.99)  GO  TO  99 

X0=XMAX=xm 
Y 0=  YHAX=Y ( 1 > 

00  98  K=2 ,N 

— IF(XtK) .LT. XG ) XG=X(K)  

IFCYCK)  .LT.  YO)  YO=Y(K) 

IFCXCK) ,GT.XMAX>  XMAX=XIK> - - 

1 F(  Y(  K)  • 5 T . YM AX)  YMAX=Y(K) 

98  CONTINUE  - - 

99  CONTINUE  < 

CX=10C./t  XMAX-XC) - ..  ... 

DY=6Q ./(YMAX-Y2) 

- - — 00  1 I*1,1C1  - - 

1 G (1 ,1 ) =AXIS 

• DO  7 1 = 1,61  - - -- 

7 G (1 , 1) =AX IS 

— DO  2-  1 = 2,101 -- 

00  2 J=2 , 61 

2 G (I , J) =3L  ANK  - 

* PLOT  POINTS 

• DO  3 1 = 1,  N - - - - - - 

* 1.  IS  AOOEO  SO  X=C  DOES  NOT  GIVE  A 0 INDEX 

IX= (X ( I ) - XO ) * DX+ 1 • 

IY=(Y(I)-Y0)*DY  + 1. 

3 G (I  X,  IY)  =POINT  - - 

DEL  X=1,/0X 

0ELY=1  ,/OY  - - 

• PRINT  6,  XO»XMAX,DELX,YO»YMAX»DELY»N 

6 - FORMAT (1H1,1GX,*XMIN,XHAX,DELX=*,3F10,3»2X,  • -• 

♦*YMIN,YMAX,OELY=*,3F10 . 3,2X,*N0.  PTS,  = *,I6> 

00  U 1 = 1,61  - - - - 

■*  STEP  NECESSARY  SO  Y-AXIS  COMES  OUT  RIGHT  SIDE  UP 
11=62-1 

L PRINT  IOC, <G(J,II), J=l,101»  ' 

•100  FORMAT  C1H  ,1GX,101A1)  

RETURN 

END  
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SU3 ROUTINE  PLCTTER(X,8F,XB,XE,X0,Y0,N0X,NNX,N0Rf1,N)  396 

OIMENSION  X(l),Oe('l)  ‘ •••■  397 

C0MM0N/YMB/VM,9  393 

IF ( NORM. £0 , ? ) GO  TO  3 399 

CALL  SLlN<CF,N0RM,F9,FS,N)  + 00 

GO  TO  2 • • * 401 

VH=1.3  402 

8=5.0  403 

CALL  BOXER<XO,YO,lC.O,lO.O,X3,XE,NC)X,NNX)  +0  + 

PRINT  10  405 

FORMAT (*  BOXER  CALLEO*)  406 

EX=(Xlll-XB)*ie.0/(Xt-X8>  4C7 

IF (EX. LT. 0.0. OR. EX. GT. 10.0)  GO  TO  25  408 

Y=YH*BF(l)+3  409 

CALL  PLOT (EX»Y»3)  410 

CALL  SYMBOL (EX, Y, .2 ,1,0 .3,-1)  411 

PRINT  9,X(1),BF(1),EX,Y  412 

GO  TO  21  - . 413 

CALL  PLOT (0,0*3)  414 

00  1 1 = 2,  N ------  415 

EX=(X(I)-X3)*10.0/(XE-XB)  416 

IFIEX.LT.G.P.OR. EX.GT.10.0)  GO  TO  1 . ( 417 

Y=YM»BF(I)*8  - 418 

- 1FCM00CI,5C) .EQ.O)  PRINT  9 , X (I  I , BF ( I) , EX , Y , I 419 

FORMAT (2(0PF12.4,1PE13.4) ,110)  420 

CALL-  PLOT  (EX,  Y»2)  - * - - 421 

CALL  SYM30L(EX,Y,. 2,1,0. 0,-1)  422 

-CONTINUE  - - - 423 

RETURN  424 

- ENO  • - - 425 
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SU9P0UTINE  30XER(X9,YC»XL,YL»X8»XE»N9X,N'IX) 

- C0MH0N/YM8/YM,P 

COMMON  /C0NST/PI,XC,ALFHA,XMINUS*X1»X2,X3,X4, XPLUS  , X JN , XN 
♦U.XM 

1 « XN,DEL  * DEL  2, A, EXPA ,XNU2 « ALFA4 , ALFA 2 » ALFA 42 
CALL  PLOT (X0,Y0,-3)  - - - 

PL0T(XL,G,2) 


°LOT (XL, YL,2) 
PLOT ( C ,YL,2) 
PLOT (0,0 ,2) 


CALL 
CALL 
CALL 
CALL 
SIGNbI.Q 

V6=C.O  - - - - - 

OY=YB*-SIGN* » 12 

- X=XL/NOX  - 

L U=  NO  X - 1 

WRITE(6,130>  X ,L  U ; — 

FORMAT (F12. 5,110) 

00  1 1=1, LU  - - 

CALL  PLOT (X,YB»3) 

CALL  PLOT ( X, OY, 2)  

X=X  +XL/NDX 

1 F(  SIG  N.L  T . 0 . 0 ) GO  TO  5 

SIGN=-SIGN 

- Y8=10 « 0 - - - -- 

GO  <0  4 

- XFI=0 « 0 — - 

WRITE (6,110) 

FORMAT  (*  X-AXIS  COMPLETEO *) 

SIGN=1.0 

- DX=XFI*SIGN*.  12 

Y=.  5 

00  2 1 = 1,19  

CALL  PLOT (XFT,Y,3) 

- CALL  PLOT (0 X, Y , 2 ) : 

Y=Y*. 5 

- IF(SIGN.LT.O.O)  GO  TO  7 

SIGN=-SIGN 

XFI  =XL  — - - - - • - 

GO  TO  6 

CONTINUE  • 

WRITE (6, 120) 

FORMAT (* 

CALL  PLOT (X0,Y0,3) 

Y=-l. 0 

X=1 .0 

00  3 1=1,9 
XI=XMXt-XB)/XL*XB 
IF(XI.LT.XC)  GO  TO  3 
SO=SQRT (XO/XI ) 

0EG=57.2958*AC0S(SQ) 

PRINT  10, X, XI, DEG, I 


Y-AXIS  COMPLETED 


427 

428 
*26 

429 
43c 

431 

432 
*33 

434 

435 

436 

437 

438 

439 
443 
*41 

442 

443 

444 

445 

446 

447 
446 
4*9 

450 

451 

452 

453 
-454 

455 

456 
45? 

458 

459 
463 
461 
*62 

463 

464 

465 

466 

467 

468 

469 

470 

471 

472 

473 

4?-, 

475 
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FORMAT  C3F12. 5, 1101 
CALL  NUM3ER(X,Y, ,1- ,O£G,0. 3,11 
X=X*1.0 
XNUM-X8 
V*-.5 

x=o  • ■ 

oo  e i*i, nnx 
XNUK=XMXE-X9)/XL*X8 
PRINT  10  ? « X ,XNUM 
FORMA T < 2F  12 • 5) 

CALL  NUM3EP  (X,Y».15, XNUM  ,0*0,1) 
X=X+XL/ tNNX-1 ) 

X=-l. 0 

Y =0  - 

YNUM=(Y-3)/YM 
00  9 1=1,11 

CALL  NUM3ER(X,Y, .15 ,YNUM, 03 . 3, 1) 
PRINT  ie,X,r,YNUH,I 
YsYU.O 

YNUM=(Y-9)/YM  - - 

l RETURN 

END  - - 


j 

i 


■ 

; 
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SUBROUTINE  SLIN(F,N0RM,F3,F3,N>  493 

DIMENSION  F ( 1 1 ••  4,99 

COMMON/YMB/YM,B  JC3 

FBsFdl  • - * 5C1 

FS*FCH  j,C2 

00  1 1 = 2, N - — - • 503 

IF(F(I) .GT.F9I  F9=F(I)  5C4 

IF(F(II.LT,FS1  FS=F  (I)  505 

CONTINUE  505 

IFINORM.NE.l)  GO  TO  25  507 

IF(ABS(FS) • GT .FB)  F B=FS  528 

00  20  1*1, N 509 

F(II=F(I)/FB  510 

8*0  .0  • — — - 511 

YM=10.0  512 

PRINT  12  ,FB,  FS  ,Y  M,  B 513 

IF(FS.LT. 0.0)  GO  TO  10  • 51V 

RETURN  515 

YM=5.0  516 

9=5.0  - ~ — • ■ 517 

PRINT  12 , FB,  FS  ,YH,  B 518 

RETURN  - < • 5ig 

IF(FS.LT.  O)  GO  TO  2 . . 520 

8=0  - 521 

U=.C0001  522 

~ IFtU.5T.FB)  GO  TO  3 - 523 

U=U*10.0  » 524 

GO  TO  4 — - 525 

YM=1C . 0/U  526 

- IF<YM*FB.LE.5.0)  YM=2.t?*YM — - 527 

PRINT  12  »FB,FS ,YM, B 523 

RETURN  . - - • - 529 

IF(FB.LT.O)  GO  TO  7 533 

•0=5.0  - - ■ — - -531 

IF(ABS<FS) .GT.F8)  FB=ABS(FS1  532 

U*.  00001  • ? -533 

IF(U.GT.FB)  GO  TO  5 534 

• U*U*10.0  --  535 

GO  TO  6 536 

YM=5.  0/U - 537 

IF(YM*FB.LE.2.5)  YM=2.0*YM  533 

PRINT  12,FB,FS,YM,B  539 

RETURN  549 

8=10. 0 - — ■■  • 54i 

FB=A9S(FS»  542 

U*. 00001  - - 543 

IF(U.GT.FB)  GO  TO  8 5**, 

•U=U*10.0  - - - - 545 

GO  TO  9 546 

YM=10  • 0/U  - •••  - - - 547 

1FC YH*FB.LE .5)  YM=2.0*YM  543 

PRINT  12,FB,FS,YM,B  549 

FORMAT  C4F12.5I  - 559 

RETURN  551 

5 NO  552 
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COMPARISON  OF  SOLAR-WIND  PARAMETERS  AND  GEOMAGNETIC  ACTIVITY  INDICES 


Philip  M.  Pazich 
William  F.  Bellew 

March  1976 

Magnetospheric  Dynamics  Branch 
Space  Physics  Division 


I.  DATA 

A comparison  of  solar-wind  parameters  and  geomagnetic,  indices  was  made 
to  determine  the  parameters  that  could  best  be  used  for  prediction  pur- 
poses.  The  solar-wind  data  consisted  of  1000  hourly  averages  of  plasma 
(v  and  p ) and  magnetic  field  data  (Bx,  Ejy,  Bz,  <Tx,  Ty,  0*z)  from  the 
Explorer  33  satellite.  The  data  were  measured  between  19  January  and 
6 April  1968.  An  earlier  study  by  Garrett  et  al.  (1974)  also  used 
data  from  the  Explorer  33  satellite  measured  during  this  period.  Approxi- 
mately half  of  the  data  used  in  the  current  study  coincided  with  the 
data  set  of  Garrett  et  al. 

Solar-wind  data  for  each  hour  were  recorded  on  magnetic  tape  along  with 
the  Aii  index  for  the  same  hour  and  the  following  two  hours  so  that  cor- 

v! 

relations  could  be  done  easily.  The  Kp  and  ap  indices  were  read  from 
an  additional  data  file  as  needed.  A given  hour  of  solar-wind  data  was 
correlated  with  the  Kp  and  ap  indices  for  both  the  three-hour  interval 
in  which  the  solar-wind  data  were  measured  and  the  following  three- 
hour  interval. 

M 

t 


II.  CORRELATIONS 


Functions  of  solar-wind  variables  were  detemined  on  the  basis  of 


physical  arguments  and  previous  studies  to  be  those  that  gave  the  high- 
est correlations  with  magnetic  activity  indices. 


The  indices  are: 

AE  0,  AE  1,  AE  2 

Kp  0,  Kp  3 
ap  0,  ap  3 


Auroral  Electrojet  index  lagged 
0,  1,  and  2 hours 

Kp  and  ap  indices  from  the  same 
3-hour  interval  as  the  solar-vdnd 
data  (Kp  j,  ap  0)  and  from  the 
following  3-hour  interval  (kp  3>  ap  3) 


The  solar-wind  variables  are: 
BZSE 
BZSM 
RHO 
v 

SIGT  or  ST 


- magnetic  field  component  Bz,  SE  coordinates  (gammas 

_ ir  it  n n sM  11  '■ 

- solar  wind  ion  density  (ir/cm3) 

- solar  wind  velocity  (km/sec) 

- (CTx  ^TyN^Jt  for  the  hourly  average 
Correlations  between  the  magnetic  indices  end  functions  determined  from 
the  solar-wind  variables  are  listed  in  Table  I.  Brackets  indicate  that 
only  negative  values  of  the  quantity  were  used;  values  > 0 were  set  equal 


to  0. 

III.  MULTIFILE  REGRESSION 

A multiple  linear  regression  was  performed  using  the  functions  and  indices 
listed  in  Table  I.  The  highest,  consistent  correlations  were  for  the 
functions  v^BZSM^  and  vSIGT  for  indices  AE  1,  Kp  0,  ap  0.  The  re- 
gression equations  for  these  three  indices  as  detemined  from  the  entire 
data  set  are  given  in  Table  II.  The  function  <BESi:^  is  the  solar  nag- 
netospheric  Bz  component  when  that  component  is  negative  and  is  set  equal 
to  zero  for  positive  values  of  the  component.  In  the  1000  hours  of  da&& 
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there  were  509  hours  with  an  average  negative  Bz  and  491  hours  with  an 
average  positive  Bz.  Additional  terms  to  the  regression  equation  did  not 
significantly  improve  the  multiple  correlation  coefficient. 

The  data  was  divided  into  two  500  hour  segments  to  determine  the  effect 
of  a specific  data  sample  on  the  regression.  The  results  of  these 
regression  analyses  are  listed  in  Table  III.  For  the  Kp  and  ap  equa- 
tions, the  solar-wind  functions  were  selected  in  opposite  order  by 
the  regression  procedure  for  the  two  data  segments.  This  indicates  that 
the  inclusion  of  the  additional  term  serves  only  to  increase  the  accuracy 
of  the  fit  for  the  specific  data. 

The  regression  analyses  for  the  A£  index  lagged  by  one  hour  always  gave 
the  highest  multiple  correlation  coefficient  and  the  computed  coeffi- 
cients in  the  regression  equation  were  the  most  consistent  between  data 
samples.  An  equation  for  Kp  and  ap  determined  by  multiple  linear  re- 
gression on  specific  data  of  v,  <BZSh>,  and^SIGT  could  be  expected  to 
result  in  correlation  coefficients  of  between  .5  and  .6  when  applied 
to  another  data  set.  An  equation  to  predict  AE  1 would  give  a some- 
what better  correlation  of  between  .6  and  .7  when  applied  to  a different 
data  set 

REFERENCE 

Garrett,  H.B.,  A.J.  Dessler,  and  T.W.  Hill,  Influence  of  solar  wind  varia- 
bility on  geomagnetic  activity,  JGK,  79,  No.  31,  November  1974,  F*  4603. 
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Table  II.  KULTIFLE  LINEAR  REGRESSION  FOR  1000  HOURS  OF  SOLAR-WIND  DATA 


R=  Multiple  Correlation  Coefficient 


Kp  0 = -.00219  v-<BZSM>  + 5.43 

= -.00216  v<BZS:,>  ♦ .00332  vSIGT  + 2.68 


R=  0.492 
R=  0.666 


ap  0 = .0117  vSIGT  +2.99  R=  0.463 

= .0115  vSIGT  - .00546  v 3ZSM  + .00319  R=  0.618 


AS  1 = -.146  v*<?ZS:>  + 114.5 

■=  -.145  v.<BZSK>+  .0863  v-SIGT  ♦ 52.3 


R=  0.699 
R«  0.732 


i 


Kp  0 = 

-.00242 

v*<BZSK> 

♦ 5.23 

R*=  0.470 

m 

-.00249 

v.<BZSM> 

♦ .00344  v-SIGT  ♦ 2.54 

ap  0 ' 

-.00572 

v-<3ZSK> 

♦ 78.5 

R«  0.450 

m 

-.00590 

V-  <3ZSH> 

+ .00817  vSIGT  + 1.47 

R=  0.621 

AE  1 = 

-.140 

v*<BZSI'>  ♦ 

102.9 

R-  0.631 

= 

-.142 

v<B2S!I>  + 

.101  vSIGT  ♦ 23.6 

R«  0.701 

b)  Second  500 

Hours 

- 

Kp  0 = 

.00510 

vSIGT  + 

3.59 

R=  0.532 

z 

.00467 

vSIGT  - 

.00183  v.<BZSM>  ♦ 2.61 

R-  0.695 

ap  0 « 

.0185 

vSIGT  + 

.0140 

Br  0.569 

- 

.0175 

vSIGT  - . 

,00459  v<PZSK>  - 2.46 

R”  0.658 

AE  1 = 

-.145 

v<3ZSM>  ♦ 

129.5 

R=  0.731 

- 

-.142 

v<BZSM>  ♦ 

.0767  vSIGT  ♦ 79.7 

R=  0.749 

■1 

t . 

L , 
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(b)  SPECTRAL  ANALYSIS 


The  second  area  of  data  handling  was  in  preparing  micropulsation  data 
for  power  spectral  analysis.  The  completed  programs  from  this  effort 

a.  take  the  data  from  the  archive  tape  and  convert  it 
into  engineering  units  and  give  a quick  plot  of  the 
data  (Program  2)  and 

b.  filter  the  data  so  lat  it  will  be  in  a form 
suitable  for  spectral  analysis  (Program  3). 

Work  is  going  on  to  produce  a contour  plotting  routine  to  display 
dynamic  spectra,  that  is,  spectral  density  contours  will  be  plotted 
on  a time  - frequency  graph. 
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PROGRAM  2 


ft 


4 


This  program  unpacks  data  from  the  archive  tape  and 
plots  it  on  the  printer  for  a quick  look. 

The  subroutine  LEVELER  suppresses  points  that  differ 
from  the  average  by  more  than  one  RMS.  This  is  to  ensure 
that  the  whole  plot  will  not  be  dominated  by  a few  large 
spurious  values. 

PROGRAM  3 

In  order  to  properly  frequency  analyze  a time  series 
it  must  first  be  filtered  to  avoid  aliasing.  We  assume 
what  we  are  analyzing  is  a continuous  process  in  time, 
however  it  is  sampled  at  discrete  intervals,  say  At. 

For  instance  let  the  power  spectrum,  P(f),  of  a 
stationary  process,  X(t),  be 


P(f)  = J*  C(  x)  e~  1U)Tdx  ( co = 2 tt  f ) 


where 


C(t)  = lim  f X(t)  X( t+x) dt 


-T 


is  the  auto  covariance.  Let  the  values  of  C(x)  be  given 
only  at  uniformly  spaced  values  of  t,  t=0,  ±At,  + 2At  ±.. 
we  can  calculate  the  aliased  power  spectrum  Pa(f)  by 

00 

Pa(f)  = J* V(x,Ax)  C(x)e"la,Tdx  (1) 

.00 

where  V(x,Ax)  is  an  infinite  Dirac  comb 
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A(x ,At) 


The  teeth  of  the  comb  are  6-functions.  Applying  the 
convolution  theorem  to  equation  (1)  we  get 

P a ( f ) = A(  f ; 1/ At)  • P(  f ) 

where  A is  also  an  infinite  Dirac  comb 
q =oo 

A( f ; 1/At)  = £ 6( f-q/At) 


Explicitly  we  have 


Pa( f ) = £ P( f-q/ At)  (2) 

q = -oo 

From  eq . ( 2 ) it  is  easy  to  see  how  high  frequency 
components  make  contributions  to  low  frequency  components 
in  the  spectrum,  see  Fig.  1 


P( f+At) 


P(f)  P( f- 1/At) 


With  this  cut  off  Pa(f)  is  equal  to  P(f)  for  f < fN=l/2Ar, 
this  frequency  is  defined  as  the  Nyquist  frequency. 

To  avoid  aliasing  the  stationary  time  series  I ( t ) must 
be  passed  through  a filter  F(t) 


4>(  t)  = F(  r)  I ( t-x)  dt 


Transforming  to  frequency  space  we  have 
♦(f)  = R(f)I(f) 

If  <|>(f)  is  to  reproduce  1(f)  exactly  for  f < f^  then  R(f), 
the  filter  response  must  be  of  the  form 


R(f) 


To  achieve  such  a filter  response  requires  an  infinitely 
long  filter  in  the  time  domain  (This  is  analogous  to  a finite 
wave  packet  in  space  requiring  an  infinite  number  of  frequency 
components).  Taking  only  a finite  amount  of  data  - letting 
the  filter  cut  off  after  some  time  leads  to  Gibbs  phenomenon, 
that  is  the  filter  response  oscillates  in  the  neighborhood  of 
the  cut  off  frequency. 

This  can  be  mitigated  by  not  requiring  a sharp  cut  off. 
This  is  affected  in  Program  3 by  replacing  the  sharp  cut  off 
by  a sine-terminator,  adding  half  a sine  wave  to  the  filter 
response 


f 

-fN  ^ 

1 fN 
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The  program  then  calculates  the  weights  Wn  = F(kA-r)  that 
make  up  the  filter.  This  filter  is  described  in  a NASA 
report  by  Behannon  and  Ness  (Design  of  Numerical  Filters 
for  Geomagnetic  Data  Analysis,  NASA  TN-D-3341). 
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Program  2 


F 


1 


PROGRAM  U^AC<(OUTPUT,TAPct,TAPfc‘.’,TAP'?,  T A r>r  i, , TflPFf  ^OUTPUT  10 

DIM:  NS!  Of.  NP  (72)  ,NU<  "NO)  .'HYd  3>  ,'J£WT  <’*>  , IDIOT  {<♦)  , IX  (16?*  11 

4) 1 1 


_ ry/My  ^ J ( i i fjJ — * — - ^ 

, ,T~~  " 

iyy  ^ * \',U  1 1 <*  1 nfjffJtl±di^^ 

^ { <t)  = j ^ tj  > * 7"> , ~~ ■ 

-ifim(ii|-t-  -■»«  ♦•«,ui  

fFlHjll . ..  J ^ 1113)  ay#  *<?_*,  , — 

>_>  * * *w  » * 31 ) . tjij777tT~^~~'-^~~ 

PwrTiTrrfi ~-~ — ‘ 1 a J > •T>ToT?TT7m777yt— r 

J ’ u,1'^U7777T~^t- _ •"'•'D, i=i5i 


i j r7ii^F>T— — 3:!II!l!^^ 


♦ US>  -1 ^h-^Ix_lVAsr  rr^T  " 

^r1 “ — LllSi~^LM^A^Ujc^  , 


SALLicuf,  _ 

cnhrTijij?  J,v-  » C;Htu~nr)p—: — - — 


=X  COORO  MATRIX 


flliuiwn 


YMAME 


MUST  MATCH  OHcNSION  OF  Y(T,J>  IN  MAIN  PROG. 

"i 


=L AST  Y FIELD  TO  P-  ^LOTED 


=6  CHARACTER  TITLE  Y-AYIS 


=MAX  VALUE  Y IF  Y>*?  PLOT  AS  Y2 


=GRAPH  TITLE  (A6  FORMAT) 


=PL0T  SIZE  CONTROL 


= 65  PRINT  COL 


(TELETYPE) 


DIMENSION  X( 1) ,Y  (IOY, 1) , A(1j1)  ,AX  IS  (6) . 000(1)  , MARK  (11) 

i?l 


wlwlm 


DATA  STAR , ILANK/1H* , 1H  / 


Hit' 


OATA  MARK  /I  Hi , 1 H2 . 1H3 , 


uutti 


IF (N3C 0 . 3T«  20)  N3CO=20 
HAI  F=  (22-NAC 


LAST =101 

MU  NO  si  j j . 1 


IFCISZ  . NE  • 1)  GO  TO  1C 
IFINNCn  .GT.  10)  MHO 0=1 T 
IHALFs(ll-N‘3CO)/2 

LASLri.  1 

HUND=6 y . 0 


H3.1H6.1H7. 1HP.1H9.1HL.1HX/ 


IF  (A  ( I Y > . -Q  . 5LAMK)  A ( I Y ) = MA=?K(  << ) 



I F (A  ( L A ST ) .TG.  DLAMO  A(UST)=ZI 

) 


irasz  .nr;,  d go  to  Z02 

. URUllb ) XiHL*J  \LL±*L=1*±£2I1 T 

GO  TO  SJ9 

-UlIILI'i..  L L'SA  ilJ  j tuX-LLL,^! 

LlfFT-i  ir.T  ♦ i 

. JLLU.  -.Ll/lLll 

if (ir.z  .1:::.  i)  go  n i j o l 

\ '■  ,i  I T L (.0  , b J a ) .<:  i A 1 1 L . 7. 1 , ( J A3  ,1 » 1 = 1 . . « ) . , 21 , ( SASh  »J_r.  1 j,. 

* zi , (hash,  r-i , ■>)  , zi,  (HASfsTri , v>  * 


61 C FORMAT (4X46. 4X51AI) 


203 


SUrtMUTINE  DECODE R (I  1,K,  IE9F) 


r ft t ft  iirixd). i*i. i2i)/i53.i6Lti7«»ies."»5-:*ior, i,5i. m 


♦ ♦ 1 0 ♦ 1 y 1 , lfc  1 » 1 f 1 ♦ 1 9 j , ;•  » £ * 1 3 c , fc 


2 1 u 6 , 9 , 5 9 , 1 J ' J t i 5t  « 1 6w  , 1 7 2 « 1 9 1 * 1 0 ♦ o f * 1 1 1 * 1 1 » f 1 * 1 1 1 » 1 2 * b Z , 1 1 


1113,  14,64,  11  4, 153,19  3,  173,1 92, 15,65, 11 5 116,17,67,1 


K = K*1 

FFAjC)  I'UFF 


mum-lam 

— p ri in  r » * sta  t i an i ♦ E2 mzzzz = " a iuzz iy,  iu  2 ud. 

JZU2i>)  -l-\ 

— IS1^=  i JA SdlF  1 1 13  ( 2 ii  l *A-  22  

IS2.=Zu*SHIFT(Ij  (23  J>  .A.  4 1,56) 


IHZ.  SHIFT  (13  (233)  .A.l<,  J!l,54) 


IH-<,Irl2a,  IHl?f  IH8,IH4,H2,TH1 


101= SHIFT ( 13 (239) . A. 2 l OB, 53) 
■ )K,n;Ti  r.->  r ■?  f i 


ID4=4*SHIFT< I M2o9). A.13 J]1t51) 


?Lo 
2 47... 
?t\ 

■ T-4  3 — 
2tJ 


) » 55 ) 


JY1=C.H]  FT  ( I.H24J  ) . A.  2 ri  33, 33) 

„JLX2 2.2+  i'UilULLUi.- IJ  . \ .h  j 13.  52)  . 
IY  ♦=-''  SO  TF  T(  I‘)  (?<>J  ) • 1 . lu  3 ? 3»  5 1 ) 
t c u r>  i 


13  (22i  > =2*  IY  if 

_L=_ix2i^ 

IFiniI).'i-.L4l  MIDsnin-itHi 


79 


5'J3-.n  JT  If  • L "V c.L c *“?  (I  <iY  i * NSS A V ) 

313 

DI  li:'  “‘SION  I*  ( i ) , Y 1 1 1 

PRINT  13c 

32  0 

POPM  \ T 1 ^ U <L  ) 

■»  ? t 

Nj  1i  = i>iS5fiV-l 

322 

« 

m-  t ^ /?!<?';  a <i 

*7  0? 

l-<  ' 

r*i 

i i n 


.iOItJO  T J S**_. 


360 
2± JL 
362 


• I 


1004*0 

10046 9 
103609 
100610 
ico«?e 
100930 
103549 
100550 
100600 
100610 
109623 
103630 
100640 


iB if i b r 

rnifzo  t 

101430  T 

IP  If a B I 

1014*0 T~ 

151700  I 

iomc — r - 

1017?0  I 

10i>30  -T~ 

lcma  i 

- isitto r~ 

101*00  I 

10  If  1 B — I - 
1018?  0 I 

toisso  — r~ 
laiff o r 

— toto«o r— 

101400  I 

- 181 41 P t~ 

ioi«?o  i 

101430 ?- 

101943  I 

—101400 T- 

1C  700  0 I 

1CZC10  T 

iczczo  r 

107033 f- 

107040  T 

— 147000 r~ 

107100  I 

107IIB  — T- 
107170  I 

18713  0 - ~T  - 
10711.0  I 

—107100 1— 

1077C0  I 

107710  “I 
107770  I 

107730  1 

107740  I 

~IO??fO 1— 

107700  I 

107310  ~t 
107370  I 

107330  f 

107340  I 

— 1373FT T~ 

10  740  3 r 

107410  ' ' T ‘ 
107470  I 

107430  I 

1C744B  I 

107470 r~ 

107*00  I 

107410  I ■ 

107*70  I 

107430  T 

107440  I 

-107440 r~ 

107*00  I 

loztio  — r~ 

107470  I 

107*30  — I 
107440  r 

107440  I 

107700  T 

107710  T 

107770  t 

107730  1 

107740  I 

— 107740 T~ 

107000  I 

10741  0 I 

107870  r 

107030  1~ 

107440  I 

172855"  T~ 

107400  I 


t 1/64 

1 T7ia 

1 1/67 

I 1,6/ 

I 1 a/ 

T i 6 T 

J 1.6/ 

I i si 

I 1 dj 


84 


lOOklO 

T ~ ~ '•  ' 

• 

I 

120.7.. 

j 10  <5 «.?  0 

I 

I 

J49,  SS 

100430 

T 

0 

1 

12*. 31 

100440 

i 

• 

I 

175.73 

V]  100450 

T 

• 

I 

12c.. la 

j 100500 

I 

• 

1 

177,-7. 

100510 

1 

• 

I 

let .50 

100570 

I 

• 

I 

177.77 

100530 

r 

• 

I 

1 it  % 1 ■ 

100*40 

i 

• 

l 

.129,1* 

" !0055(T 

i 

• 

I 

128. e? 

1006C0 

i 

• 

JL 

174.-1 

100510 

i 

• 

I 

176.45 

100570 

i 

0 

_Jt 

J29.0L_ 

100530 

i 

0 

I 

100649 

i 

0 

I 

no. on 

-looeTo- 

i 

0 

I 

1*>U.  . * 

100700 

i 

0 

r 

J 71, 1.0  _ 

' 100710 

I 

0 

i 

131.7- 

100T70 

i 

9 

i 

i J 1 * to 

100730 

I 

9 

i 

131.44 

100740 

T 

0 

r 

1:2. On 

"10  0790“ 

I 

« 

i 

112.13 

100500 

I 

* 

i 

112. 36 

iooeio 

Y" 

• 

i 

133.6/ 

100070 

i 

• 

134. . 4 

100*30 

I 

• 

U-.ol 

100*40 

i 

• 

i 

ns. a 

II 705  IT 

i 

• 

i 

132. / 1 

100500 

I 

• 

i 

137.96 

100510 

T 

0 

i 

133. .3 

100970 

i 

9 

i 

133.45 

100930 

I 

9 

i 

100940 

i 

9 

i 

lll.lo 

"1009*  O' 

i 

”1 

130.54 

101C00 

i 

0 

i 

loi. 1 o 

101010 

i 

0 

i 

132.39 

101070 

i 

0 

i 

U2.  it* 

101030 

T“ 

0 

i 

ioi.ui 

1010*0 

r 

9 

i 

1 0 . o o 

— toiofir 

i ■ "■ 

i 

1-9.59 

101100 

r 

0 

i 

179.3> 

101110 

T 

i 

1 2 o • i‘c 

101170 

r 

0 

i 

126. Jt 

101130 

r ■ 

i 

1*  6.  >1 

101140 

t 

0 

i 

179.71 

102010 

107b?0 

■ 

10790? 

101MI 

ie?«70 

»?!» 

io’s*e 
107*00 
10»19 
107*20 
to?*30 
1101(0 
197999 
197799 
19771# 
1977?# 
19773# 
10771,0 
‘ 107»90' 
197000 
107010 
107*70 
107*30 
1079b# 
~!07«e 
107*00 


r 


112. 

lib. 

lib. 

113. 


61 

t« 

71 

rs 


113. 
113. 
11  J. 

lib. 

11c.. 

lib. 


1 1 o . 
117. 
11*. 
11*. 


Ho 
It  ... 
117. 
_IU. 
no. 
112. 
ni. 


c i 
9*_ 
• 0 

to 

0 .< 

t.b 

b«T 

*u 

«* 

13 

91 


i 1 

9<. 

S9 

3i 

b0 

3b 

9v 

00 

lb 

7b 


j 


j 


11b 

L 117. 

no. 

lib. 

112. 

Ali 


I 


113. 

113. 


01 

7o 


C7 


iokoo  t 

101(10  T 

101(70  I 

101623' T~ 

101640  t 

iokoo r- 

101700  i 

loiTio  — r- 

101770  T 

1017*9  * I 

ioi7<.o  r 

1017*0 T- 

101(00  I 

101*10  I 

101070  T 

101830  ~I~ 

101*40  I 

--mm r- 

loioeo  r 

101910  1 

101970  I 

10197?  T 

10  1940  I 

- 1019*9 r- 

107000  : 

107010  T 

1C7070  r 

1C2C29 I ' 

107040  1 

— .--O**— 7“ 

107100  I 

107110  T 

107170  I 

107120  — r~ 
107140  I 

—TO  7190 r- 

107100  I 

107710  — r 

10777C  I 

107770  T 

197740  I 

—1077*0 1- 

107700  I 

107310  -1 

107370  T 

K7330  — T~ 
10  7340  I 

— 107790 1- 

107409  I 

107419  T 

107470  T 

107430  T 

107440  I 

-'107455 r~ 

10  750  C I 

107519  — T- 
107570  I 

107530  T~ 

102549  1 

—107950 T 

102(00  I 

10261 0 T 

107(70  I 

107670  I 

1*7640  1 

107**0  1 

iovro  i 

107719  I 

10777*  I 

197?30  T 

117*40  I 

1077*9  T 

19  7*0*  I 

102(19  T 

107*7?  T 

197*19  * 

1*3*49  t 

197*59  T 

10  790  0 I 


__9» 1 HZ  • 

9-  l J . " 

i-  X’j2  • 
>«1  i/2  • 

-rrt  rr 

j • 1 jZ  • 

9- 1 >4. 
~9-l3d.  ‘ 

>- 1 j‘*  . 

$ *li»o. 

Z.'  1 "h* 

is. 

9-139. 
”158 
- i;<?* 

3-  I J M . 
9-13S. 
S-loOl 
j ioL. 

~T~urr~ 

_5-;  lO  v •_ 
”9  loo . 

t-  luO. 

“ >“loi. 

9-  Id 2 
3 * io  3 • 

_9-low._ 

J lo  > • 

9 - loo  * 

3- 16c  • 

3 log . 

9 loi# 

_3  16; 

3**  lfcc  * 

1 93  • 

3- lot . 

3 lo  7 . 

3 - loo. 

3 . loO  . 

$ io7. 

_3  ;1*>  - * 

~ 3 •*  161  . 

S-lo>. 

9-  loc  • 
_j>.  Ion 
■J  loo  • 
_S^loS. 

5 • lo?* 

3 130 
3 1 O 7 • 

5*l°v . 

5 • 16 9 • 

3- lot . 

3 • loo • 

3 lod  . 

S* lo8  / 

S • lcS. 
_j  i/?. 
3-173. 
d- 173. 

> 1.3. 

3-17*. 

-JLrl^i'f  . 
3 1/1?. 

3-1/3. 

-jL'ir**— 
3 1/3 

3-1/0% 

3-1/7. 

.JuWju.. 


IIWIM  ////  FNO  or  LIST  ✓✓// 

////  isc  or  list  an 
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Program  3 


PROGRAM  PEAj  T°  I I.N:‘UT  • !'U  I VJ  r , T AP-'i  , J JW}.'?  ) 
cc  j s -•  a :i  ?-ai;  taje:  ■.;••••  it  ~.i  ta p-'i  cccceg 

DI  TENSION  WT  17  i)  1 1 NX  ( J 7 T j ) ,0UTX{-  . j) 

. common  mt<  is  n ,i>  ;a g , jx ( : ^ ) 

REAL  I NX 

DATA  KEOF,LI,):C,f.PPM,i..AGT  MC/?*",i;?4S,1/ 

DATA  L I N X , LI  N f , L I‘ I Z t lCU  f X , L OUT  Y»  LOUT  Z i < ,*>'Y  , KPZ , NT  0/10* 


DATA  LO,  LFIL  ,IrJt:C,  XMAX,XMI  I, YMAX , YMIN, Z MA X, ZM I N/45 3 ,49,  It 

» l . , 1 . 'J , , •♦ » . . / • 

FORMAT (iHl) 

I?EA  )=  j 

IP=1 

MV  A.<  = 1 

NPTS<IP=i 

LNOR<A  = ?|3P3»LFIL-l 

REWIND  1 

REWIND  2 

CALL  REAL  W (W ,WT, LFIL  , IP, i. ) 

1=3 

I = m 

F E AD  ( 2 ) I M X ( I ) • 

j.  F (E  OF  ( 2 i J c , 1 

PRINT*-  , **  EOF  (2  ) ** 

HP-»R=I-1 

NPSAV  = tlPPP 

PRINT*,"  NPPF,NP3AV=M,NPPR,NPSAV 

■ PRINT  5,  ( I N < (I) , 1=1,  (nPR ) 

IREA J=IRFADf 1 

CALL  *•  I L T E 3 ( I ’ X i f ? P-: , OU  T X , L OU 1 X » X M A X , > M T N , K PX  , LFIL  , IDFC  ,L 
+ASTREG) 

CALL  KILTER! IMX , NPP j , OUT X , L OUT X, L F I L , IQ r0 ) 

PRINT  lit 

PRINT  LOUTX=",  LOUTX 

PRINT  5,(1  NX (I ) ,1=1 , J PS  A V ) 

. PRINT  111 

PRINT  (OUT  < C I ) ,1  = 1 ,NPaAV> 

FORMAT  (1CF1-;  ,3) 

CALL  CKRATCM  (OUTX,  LO'JTX) 

FORMAT  ( //) 

IF  d Vi  AU.EO.  LASTREC)  STOP  107 

NP^Riy  ■ 4ft 

GO  TO  44 

PRINT  14'0__  ___  __  _ 

F 0 <M  A I ( v L'UF  ENGUUN 1 ; R'.  0 ON  lNPU  l TAR" ' ) 

STOP  IjIC 


CONI  I'l'J'I 

I'KlMt  1K1.(H‘I(L)  »L=1»132) 


F U !\  1 (2X,  li 241) 

co'UTMur 


Pn  IN [ ijl.2 

JljJ.iLXlUJi 

Pr  INT  1 ,'C  jf  n;<YI  ,8IG/,NTNYI  ,SMALLY 

♦Ut  in  * , I 


SUJ^  TUTINE  IFLTt.K  ( I 1 , 1 1 N , 00  T , L OUT  , X M AX , Y M IN  ,K  R . LF  , in  t C , L A 


JE-(LIN-LF)/ IJtC+1 


DO  5>5  J=l,Jz 


IF  (I  W • t'4»M  Ii)  FID  NOFILPT  = I 
. r.  :T ,12)  30  TO  4 3 


GO  TO  4 2 

IF  (IW.UE.LF)  STOP  1) 

CO  (T INUE 

FT  1=  ^T(TW)  *C  1*1  U)  » +FTO 


LOUT  = L0 JT ♦ 1 

IF  (L  Oti  T . G T , . . ' ) STOP  11 


OUT (LOUT ) =FTO 


GO  TO  i 4 

• ■ ■ : 


WRIT  F ( i ) LOUT, (OUT <<> ,<=i, LOUT) 

PRI\’  T 77  , I RE  AO  * K K , LO’J  T • '<  NA  X , X MIN , (OUT(L  ) ,L  = 1,  IT) 


IF('10D(<4,6)  .-•}.?)  CALL  S<R  ATCH  ( CU.T , 5.  " ) 
tRINT  *,  11(12-2)  ,IN(It>l)  ,IN(IE>  


XHAX  = -03  v'  -i  4'i  J43-j's)(4. 

LLIi'ir  FO  = NL_lNi±J*  1 jEC 


PRIMT * , “ LIME  COUNT  AT  CNO  OF  FILTERED  TATA 


I - ♦ I d t 


9 


AD-A035  941  EMMAN'JEL  COLL  BOSTON  MASS  F/G  17/6 

INVESTIGATION  OF  MICROPULSATION  ACTIVITY.  1.  MAGAF  SYSTEM  AOOIT— ETCIU) 
AUG  76  W F BELLEW*  C J CANTOR*  M P HAGAN  F19628-76-C-0013 

UNCLASSIFIED  SCIENTIFIC-1  AFGL-TR-76-0244  NL 

END 

DATE 
FILMED 

3—7: 


. 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  8UACAU  Of  STANOANDS  - lt«S  “ A 


N=LF IL/2 
PLP=1  ./FLOAT  (N> 


H=i./N 
N°l=  N + 1 


TWODT=2.*OT 
PI=?  .I4i59265~5393 


W(i)  =PLp 

19 

P3UM= l . 

is: 

F=C3$(H*K*PI  )/  (1  *!t*H*<  + K) 

GO  TO  * 


F = .7  i?39<?i634 

print*,"  i >i  suiv  ^eaqw 


<*pd/  (pi+  <> 

0)  PRINT  * , F , S 


W ('<♦!  ) = F*S 

psu.v=rsuN+w(,<+-i) 


CO  NT  irj'JE 

SU  10-LS  = W ( 1)  + >.*P3LM 


OELT  A = 1 . -jUiOFLS 

OF LT  A D = ")  E L T A / F L 0 A T <?*N+1> 


P»i»T*,”  SUM  OF  UNCOWcCTFJ  FILTER  WEIGHTS 


19-i 

K 

2c: 

UMOFLS 


DO  11  I 


SUMWT5=3U1WT3+W(I> 
WT  (I  «■•!)  aw  < I) 


IC=ND1* t-I 
VtT(n>=H(I) 


PRINT*.  *•  SOM  OF  CORRECTED  FILTER  W£IGHTS=",  2.  * SUMWTS- W<1 ) 


CALL  3K  <A  TCil  ( NT  » l'*  N*  1 ) 
PRINT  13i,  UPl.WT 


FORMAT  (lHi.Ia,*  L3J-PASS  FIL.  WE  IGHV3, ''  I , W 1,  . . . , VN*  , / ( IX 
♦1CE1 3. f ) ) 


no  5 ] 1 = 1, 13  ) 

SPaFLOAT  <1-1 )/lc9.» 
F*.Eu  = SP/T  ROOT 


DO  K = 1 » N 

K. 

SUMC  = i>UMC ♦ WT  (NPl+<O*C0S(PI*FK*SP) 

C oraiiiiJE 

R 

p 


0 Tfl  IJl) 


CALL  3KRA  T CH (F REORL i , 1 3 j ) 


F T J=  - T » W f (I  l)*u:i 
LOUT  =LDJT  + 1 
OUMLOUT  )=FT I 


IFiX'IrtX.LT.FT:;)  X'lAXsFTD 
IF  (X  fl I N « C f « F TrJ ) X ‘1 1*1  = FT 0 


ID=I  2+ IOcC 
IF-ir»I3F.C 


P.-5IVT  k , tounit , 1 = 1,  LvUT) 

FOriATlUFL  ,?t ■ 

F:FTJF.‘J 


25  i 
251 


?5  7 


FNO 


3.  PROGRAMS 


(a)  In  Progress 

Work  is  going  on  to  recognize  spurious  points  in  the 
data  and  eliminate  them  before  the  data  are  analyzed. 

(b)  Programs  Planned 

The  next  major  area  of  concern  will  be  to  develop 
plotting  routines  for  the  Varian  computer.  After  that 
we  will  begin  to  get  magnetic  indices  from  the  network 
data. 


